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TRANSLATOR'S PREFACE 


This book is a translation of L.A. Chernov's "PacnpocTpaHeHue BOuJH 
B Cpere Co Cnyuaiinrmn HeorHoponHOctTaAMH" , literally "Propagation of 
Waves in a Medium with Random Inhomogeneities". It is hoped that Chernov's 
book (together with the translation of V.I. Tatarski's book on wave propaga- 
tion in a turbulent medium, also in this series) will furnish a comprehensive 
ani authoritative survey of the present state of research in the field of wave 
propagation in random media, with special emphasis on important Russian con- 
tributions. 

The translation departs from the original in one small particular: the 
Russian words for amplitude and for level (amplitude measured in logarithmic 
units) have both been reniered as amplitude, since in every case the context 
precludes confusion. 

There is no ideal system for transliterating the Cyrillic alphabet into 
the Roman alphabet. Of the various systems which exist, I prefer and have 


used that of Prof. E.J. Simmons of Columbia University. 


AUTHOR'S PREFACE 


This monograph contains a systematic treatment of the theory of 
wave propagation in a medium with random inhomogeneities. Practical 
problems of acoustics, optics and radiophysics have stimulated and made 
urgent the need for such a treatment. 

In Part I we study the problem of wave propagation using the ray 
approximation. In Part II we deal with the diffraction theory of wave 
propagation. In Part III we examine the question of how fluctuations 
in the incident wave affect the diffraction image formed by a focusing 
system; this question is of considerable interest in hydroacoustics and 
astronomical optics. Some theoretical deductions are compared with 
experimental data. (The most involved calculations are relegated to 
appendices.) 

The author takes this opportunity to express his appreciation to 
L. M. Brekhovskikh, Corresponding Member of the Academy of Sciences of 
the U.S.S-R., for his constant interest in this book and for the valu- 
able advice he gave while it was being written. The author recalls 
with gratitude the late Professor G. S. Gorelik; the latter undertook 
to familiarize himself with the manuscript and made numerous helpful 


critical remarks which led to its improvement. 
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INTRODUCTION 


Inhomogeneity is a characteristic property of every real medium. Two types of inhomo- 
geneities are observed, regular and random. Regular inhomogeneities are produced by spatial 
variation of the mean characteristics of a medium, while random inhomogeneities are produced 
by deviations from the mean values. 

Both regular ani random inhomogeneities can influence (although in different ways) 
wave propagation in a medium. For example, such phenomena as refraction and duct propagation 
are caused by the influence of regular inhomogeneities, while the scattering of waves is 
caused by random inhomogeneities. The scattered waves are superimposed on the incident wave 
and lead to amplitude and phase fluctuations of the resultant field. This explains the loud- 
ness fluctuations of sound propagated in the atmosphere and the sea, the twinkling of stars, 
and (in part) the fading of radio Signals. Regular changes in a real medium occur in some 
preferred direction (depth in the sea, height in the atmosphere). Therefore, as a first 
approximation, a real medium can be regarded as having a layered structure. The theory of 
wave propagation in layered media has been developed in the last ten to fifteen years in the 
extensive investigations of Brekhovskikh 32S: 

In spite of the fact that the twinkling of stars is a problem several thousand years 
old (it was of interest to Aristotle) the development of the theory of wave propagation in a 
medium with random inhomogeneities has likewise occurred mainly in the last decade and is due 
to a large extent to the practical requirements of hydroacoustics and radiophysics. In the 
early theoretical work of Krasilnikov [40,54] and Bergmann [14] the problem of wave propagation 
in a medium with random inhomogeneities was studied in the approximation of geometrical optics. 
Applying the method of small perturbations to the eikonal equation, these authors obtained 
formulas for the mean square fluctuations of amplitude and phase as functions of distance. 

An unimportant difference in the formulas for the amplitude fluctuations results from the fact 


that Bergmann considered a spherical wave (point source), while Krasilnikov carried out the 


calculation for a plane wave. 


The author [46] suggested another method of solving the problem of ray propagation 
in a medium with random inhomogeneities, a method based on application of the Einstein- 
Fokker-Kolmogorov equation. The propagation of a ray is regarded as a random process without 
after-effect (continuous Markov process), in which the role of time is played by the path 
length traversed by the ray. The angular distribution function of the rays obeys the 
Binstein-Fokker-Kolmogorov equation. This method permits one to determine not only the angu- 
lar but also the linear displacement of the ray from its initial direction. Moreover, this 
can be done without requiring the displacement to be small, which is an advantage of this 
method over the method of small perturbations. Kharanen [5] has also used the Einstein- 
Fokker-Kolmogorov method in treating the problem of ray propagation in a medium with random 
inhomogeneities. However, his results are valid only for small angular and linear deflections 
of the ray, and in this limiting case can be obtained from the more general formulas of our 
work. 

The problem of the intensity of radiation scattered by random inhomogeneities was 
examined in early papers using the wave approach. Pekeris [35] derived a general formula for 
wave scattering, assuming that as a result of fluctuations the refractive index deviates only 
Slightly from some mean value. In acoustics it is usually necessary to take into account 
density fluctuations as well as refractive index fluctuations. In this book, we derive a 
more general formula which takes this factor into account. We also clarify the conditions 
under which, in the acoustic scattering problem, density fluctuations can be neglected in 
comparison with refractive index fluctuations, i.e., we show the range of applicability of 
Pekeris' formula in acoustics. 

In recent years, the problem of amplitude and phase fluctuations when waves are propa- 
gated in a medium with random inhomogeneities has been studied by many investigators [12, 18, 
20% 21. 2a 25% 32] from a diffraction theory point of view. In particular, it has been 
possible to ascertain the limits of applicability of ray theory for a medium with random 
inhomogeneities [12, 31, 33]. All investigators have applied the method of small perturbations 
to the wave equation; only Obukhov [12] has made use of Rytov's method ERuP which was used by 


Rytov in studying the diffraction of light by an ultrasonic grating. The fluctuations calcu- 


lated by this method are not restricted by the condition that they be small; this is the advan- 
tage of Rytov's method over the method of small perturbations. In the present book, the 

author also uses Rytov's method and, in doing so, generalizes the formulas for mean square 
amplitude and phase fluctuations obtained by Obukhov. Unlike Obukhov, who started with a 
GausSian correlation coefficient for the refractive index, the author has succeeded in deri- 
ving general formulas which are not a consequence of any special form of the correlation 
coefficient. 

Mean square amplitude and phase fluctuations are not enough to completely characterize 
the statistical properties of a wave field. The statistical properties of the fluctuations 
of a wave field can be more completely characterized by using correlation functions. However, 
the problem of correlation functions for fluctuations of the values of the basic characteris- 
tics of a wave field has not beer investigated at all from a wave point of view. In this 
book, the author tries to fill in this important gap in the theory of wave propagation in a 


medium with random inhomogeneities. With this aim, we calculate: (1) the cross correlation 


coefficient of amplitude and phase at one receiving point, (2) the autocorrelation coefficient 
for amplitude (or phase) at different receiving points and (43) the time autocorrelation co- 
efficient for amplitude (or phase). Nobody has heretofore inquired as to the nature of the 
cross correlation of fluctuations of amplitude and phase at a receiving point. In this book, 
we answer this question and show that the correlation between amplitude and phase fluctuations, 


which exists at small distances, vanishes at large distances. 


The question of the autocorrelation of amplitude (or phase) fluctuations at different 
receiving points has been studied by many authors [14,6] but only in the ray approximation 
and only for the case where both receivers lie in a plane perpendicular to the direction of 
propagation of the wave (transverse autocorrelation). In this book, the problem of longitu- 
dinal and transverse autocorrelation is studied from a wave point of view. It is shown that 
at all distances the transverse autocorrelation between amplitude (or phase) fluctuations 
extends to approximately the same separation as the correlation between the random inhomo- 
geneities of the medium itself. Moreover, it is shown that the longitudinal correlation 
extends over much greater separations than the transverse correlation. If the separation be- 


tween receivers does not exceed the distance within which a ray treatment is suitable, the 


fluctuations of amplitude (or phase) are practically completely correlated. 

Since evidently the time variation of the inhomogeneities is due mostly to their motion, 
in this book we also examine the question of how the motion of the medium as a whole affects 
the time correlation properties of the wave field at a receiving point. The correlation the- 
ory developed in this book leads to the conclusion that the time autocorrelation function is 
practically independent of distance and depends mainly on the speed with which the observer 
moves. Experiments performed on ships moving in the ocean confirm this deduction [23] . The 
theoretically calculated mean autocorrelation curve is in good quantitative agreement with the 
experimental data. 

Fluctuations in a wave incident on a focusing system are accompanied by fluctuations 
of the diffraction image. In the diffraction image we observe not only deviations of the 
intensity from the mean distribution, but also the mean distribution itself depends in an 
essential way on fluctuations in the incident wave. In this regard two problems arise in the 
theory of focusing systems: cz) investigation of the dependence of the mean distribution in 
the diffraction image on fluctuations in the incident wave and (2) investigation of the dis- 
tribution of fluctuations in the diffraction image. As far as we know, neither of these 
questions has been studied previously. Both questions are studied in this book. 

In the presence of strong fluctuations the mean distribution differs greatly from the 
distribution in the absence of fluctuations, if the dimensions of the diaphragm are larger 
than or of the same order as the correlation distance of the incident wave. The distribution 
curve falls off monotonically as we go away from the focus and, in the absence of fluctuations, 
forms a system of decreasing maxima. The distribution of fluctuations in the diffraction image 
coincides with the mean distribution in two limiting cases: (1) the fluctuations are small and 
the dimensions of the diaphragm are small compared to the correlation distance and (2) the 
fluctuations are large. If the fluctuations are small but the dimensions of the diaphragm are 
larger than or of the same order as the correlation distance, the distribution of the fluctu- 
ations does not coincide with the mean distribution. For the relative fluctuations at the 
focus the theory gives a result which coincides with that of Krasilnikov and Tatarski [10] . 

All the problems studied in this book are based on the use of the scalar wave equation, 
as in acoustics, but in passing we indicate the polarization corrections which must be intro- 


duced in going over to the case of electromagnetic waves. 


Part I 


RAY THEORY 


Chapter I 


STATISTICAL CHARACTERISTICS OF THE MEDIUM 


In a medium with random inhomogeneities the index of refraction is a random function of 
the coordinates and the time. The changes in refractive index in the sea and in the atmos- 
phere are usually caused by temperature fluctuations. Jn addition, salinity fluctuations 
in the sea and humidity fluctuations in the atmosphere can play a role. Apparently, these 
fluctuations are small in the majority of cases. Liebermann [4] has made an experimental 
study of temperature inhomogeneities in the ocean using a fast-acting thermometer mounted on 


a submarine. The curve of temperature fluctuations obtained by him is shown in Fig. 1. 
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Fig. 1 Record of temperature fluctuations in the ocean. 
(After Liebermann) 

The choice of a submarine instead of a surface vessel was dictated by the need to avoid 
pitching, which introduces additional errors into the measurements. The submarine moved 
at depths of 40 to 60 meters. At these depths the mean temperature fluctuations amounted to 
0.04° Cc, with a ~60 cm as the mean size of the inhomogeneities. Such small temperature 
fluctuations correspond to small fractional changes in the index of refraction; the mean 
square fluctuation of the acoustic index of refraction is equal to ta = 5X 107 
( vis, x 10°). Even if the temperature fluctuations amounted to a few degrees, the 


fractional changes in refractive index would still not exceed 0.01. Apparently, the influ- 


ence of salinity fluctuations is even less significant. 


1. The Correlation Function. We shall assume that the fluctuations in refractive index 
represent a random process in space and time, described by the random function of coordinates 
and time u(x,y,z,t). Regarding this random process as stationary in time, we shall character- 


ize it by the correlation function 
Ny 5 7 W(X Vp 921 t)H(%55F 53 25st) , (1) 


where the overbar designates averaging with respect to the time t, or, because of the ergodic 
hypothesis, averaging over the ensemble of realizations corresponding to the different 


possible states of the medium. Understanding the overbar to mean time averaging, we write 


+T 


a 


Bel gos 
12 = lim on U(x Y > Z45t)U(%55Y55 Zo yt dt ° (2) 
[00 “1 


For a spatially homogeneous process the correlation function depends only on the coordinate 


differences x = Xo ~ x,» y= Vo = Vp De Bee = Za Lees 


Ny 5 = Nj o(x,y,2) : (3) 


O the function No, achieves its maximum Noy 
tuation of refractive index No, = is The correlation coefficient N is defined as the ratio 


For x =ye=z equal to the mean square fluc- 


of the correlation function Ny 5 to the mean square fluctuation = i.e. 


N 
12 
oe (4) 
oC 
UL 
so that 
oO 
Nino =H N. (5) 


As the distance between the points is increased, the correlation coefficient decreases from 
its maximum value of unity and becomes small compared to unity at a distance called the 


correlation distance, i.e., the statistical dependence between the fluctuations disappears. 
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If the properties of the medium are not Spatially homogeneous, then the correlation function 
will depend not only on the coordinate differences but also on the coordinates themselves. 
However, in what follows we shall consider only the statistically homogeneous case. 

In addition to the three dimensional spatial correlation function (1), we introduce the 


four dimensional space-time correlation function 


Noo 7 u(x) 54925 9b) ulXS 5953259 t5) by (6) 


setting t, = ty + t and taking the overbar to mean time averaging, we write 


+T 
: 1 
Noo T: = oT ¥ M2) Vy 9 Zy ty )M (X51 V os Zoe ty + tat, - (7) 
=f 


For a stationary and spatially homogeneous random process, the four dimensional correlation 


function depends only on the coordinate differences and the time difference t = th . ty: eee 
Noo = N, o(x5y,2, +t). (8) 


The corresponding correlation coefficient is defined by Eq. (4). As t grows the correlation 
coefficient decreases and becomes small compared to unity for a time t ~ T called the corre- 
lation time. The three dimensional correlation function is the special case of the four 
dimensional one for ty = th: 

In a brief note Fine [24] considers the problem of averaging not only in time but also in 
space when defining the correlation coefficient for refractive index fluctuations. However, 


the equivalence of these two averaging operations for a stationary and spatially homogeneous 


process cannot be doubted. 


e. Determination of the Form of the Correlation Function. It seems that only one case 
can be given where the correlation function can be determined theoretically. This is the 


ease of homogeneous isotropic turbulence. The field of temperature fluctuations caused by 


turbulence was investigated by Obukhov (ai, who found for the mean square temperature differ- 


ence at two points a law similar to the "two-thirds law". However, even if (under familiar 
conditions) the state of the atmosphere can be described in a satisfactory way by the theory 
of homogeneous isotropic turbulence, the question of whether this theory is applicable to the 
ocean has yet to be settled. In hydroacoustical investigations it is appropriate to start 
with correlation functions found empirically. Moreover, in atmospheric acoustics (as well as 
in optics and radiophysics) it is also appropriate to begin this way if the conditions for 
the applicability of the theory of homogeneous isotropic turbulence are violated for one 
reason or another. 

The correlation coefficient for temperature fluctuations in the ocean was determined 
experimentally by Liebermann [4], who used a correlator which first multiplied together two 
temperature curves T(x, ) and T(x, + x) and then averaged them over the entire record. [In 
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Fig. 2 Correlation coefficient of temperature 
fluctuations in the ocean. 
(After Liebermann) 


the points indicate the results found by Liebermann for different distances x. The depen- 
dence of the correlation coefficient on the distance is satisfactorily approximated by the 


function 
N(x) = eo X1/8, (9) 


if we set the correlation distance a = 60 cm. Fora statistically Lsotropic medium we can 


write 


N(r) = ela (10) 


instead of (9), Where r = c= + y" + or 


A peculiarity of the correlation function e” Ix |/2 should be pointed out, namely, its 
derivative at x = 0 differs from zero. This is possible only in the case where the refractive 
index fluctuation u(x) is a discontinuous function. Actually, we assume the opposite, i.e., 
that the refractive index fluctuation u(x) is a continuous function of x. Differentiating 


the product ux, u(x, + x) with respect to x, we obtain 


= [ux, u(x, + )] = ule,) Sule, + x) - (11) 
Since 
Bey +) = a (x, +x), 


in the right hand side of (11) we can differentiate with respect to x, instead of with respect 


tO 7x 


S fa(x, u(x, + x)] = w(x,) om uc, + x) . (12) 


Averaging (12) with respect to the coordinate x, we obtain 


1 
an, (x) —————— 
= — = u(x,) ae HO +) (13) 
where 
+X 
ar ay ae dl. 
Nyp(x) = Ga )aGq +x) = Lim ge fwGn ubq + x)eR, 
=x 


We now find the value of the derivative at x = 0: 


— ae) x 
Pat =a) = *) ose a ; Ge. te 
xX=0 


dx 


Since u(x) is everywhere bounded, the limit of the last expression is zero, i.e. 


an. (x) 
Poe a a) 


The correlation coefficient -" |X| /2 does not satisfy this condition. Consequently, the refrac- 
tive index fluctuation u(x) is a discontinuous function for this case. Actually, the temper- 
ature discontinuities. are smoothed out as a result of heat conductivity, so that there is not 


- |x| /a 


a full correspondence between the correlation coefficient e and the actual conditions. 
To achieve such a correspondence the function e” |x |/2 must be modified in such a way that the 


new function differs from the function oll /a only in the immediate neighborhood of zero and 


has a vanishing derivative at x = 0. For example, the function 


piaa, g*|xifa, 2. 22 a)x| 


N(x) 2 aa - 1 aa - 1 


satisfies these conditions. The smallest value of a corresponding to Liebermann's experi- 


mental data is 0.5 on 


The experimental data are also described in a satisfactory way by a function of the form 
~x°/a” 


Which is especially convenient for theoretical investigations, and has a vanishing derivative 
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at x = Q. The corresponding correlation coefficient for a statistically isotropic medium has 
the form 

2) 2 
at’ /a 


N(r) = (15) 


In the following we shall use the correlation coefficients (15) and (10) in studying 
special cases. In this regard it should not be forgotten that the correlation coefficient (10) 


corresponds to discontinuous changes in the refractive index fluctuations. 


Chapter II 


RAY STATISTICS 


We turn now to a consideration of ray propagation in a medium with random inhomogeneities, 
assuming that a, the scale of the inhomogeneities, is large compared to the wavelength X. In 
hydroacoustics this condition is often satisfied for ultrasonic waves; in the atmosphere the 
condition is met for light waves, since the inner dimension of the turbulent fluctuations in 
the atmosphere is of the order of 1 cm [49]. It should be noted that the condition that the 
wavelength be small compared with the scale of the inhomogeneities is only a necessary condi- 
tion for the geometrical approximation to be suitable, but not a sufficient condition. If this 
condition is met, then the ray theory can be used in regions of linear dimension L, where L 
satisfies the condition Vil << ae This condition has a simple physical meaning: the size of 
the first Fresnel zone for the distance in question must be small compared to the scale of the 
inhomogeneities. At larger distances which do not satisfy this condition the ray approxima- 
tion cannot be used, and in this case diffraction theory is necessary. 

In this chapter we shall restrict ourselves to the ray model and we shall assume both of 
the conditions ’4 <a and veue << a. The necessity of the second condition will be rigorously 
justified later (Section 21). Moreover, we shall assume that the transit time of the ray is 


small compared to the characteristic scale of changes of the inhomogeneities in time. 


3- The Ray Equation. The ray equation in the form most convenient for the considerations 


to follow can be obtained from Fermat's principle 


B 
OO 2. se 
[ ~ = min. (16) 
A 


Introducing the refractive index n = ce /e; Eq. (16) can be written in the form 


B 
J aly, 2)a0 = min. (17) 


A 
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We shall assume that the ray trajectories belong to the family of curves expressed by the 
equations x = x(u), y = y(u), z = z(u) and passing through given points A and B. The para- 
meter u is supposed to be chosen so that it takes fixed values Uy and Us at the points A and 
B. For example, any of the three Cartesian coordinates satisfies this condition. The arc 
length 6 of the curve does not satisfy this condition, since it changes in going from one 


curve to another. Since do = xt Hy ohne Zr du (where the prime denotes differentiation with 


respect to u), the Fermat principle can be written as follows: 


“2 
[ nC Kyle) V8 4 wine ae du = min. (18) 


a 


By introducing the parameter u the problem is reduced to an ordinary variational problem. 


Designating 
+ a C 2 
POG yee ay ee) = nx y2) xi ty’ + 2” , (19) 


we write the Euler equations of the variational problem as 


e| 2 
“Hy 
Ne 

I 
yy 
i 
Oo 


Because of (19) they take the form 
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and returning to the variable o, we rewrite Eq. (20) as 


a(ns, ) on F 
do ~ 30x ? 

d(ns. ) - on _ ; 
do oy “ 

d(ns_ ) ; dn : , 
do OZ” 


Instead of the three ray equations (21) we can write one vector equation 


d ns) 


| i ie 
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If the refractive index is given as a function of the coordinates, then Eq. (22), together 


with the equation s = dr/do, allows us to find the equation of the trajectories r = r(d) 


with given initial conditions. 


4. The Ray Diffusion Coefficient. We shall assume that the index of refraction deviates 


only slightly from a mean value equal to unity, i.e. 
nC Ky 32) =1+ xyz) ) | | <l. (23) 
Then Eq. (22) can be rewritten as 


d(ns) 
dg 


-~ W=0. (24) 


We now calculate the mean square deviation ee of the ray from its initial direction after 
going a distance Ad. We choose the path Ad so that it is large compared to the correlation 
distance of the refractive index, but so that the deviation of the ray along the path is still 


small. Integrating (24) along the path Ao, we obtain 


n'g' - ns = f VW do . (25) 


It is physically clear that the angle of deviation of the ray from its initial direction is 

determined by the change of the refractive index along the whole path Ao, and that the longer 
the path Ao compared to the correlation distance, the less the deviation depends on the values 
nandn' of the refractive index at the ends of the path. Therefore setting the random values 


of the refractive index at the ends of the path Ao equal to the mean value (unity) we obtain 


AS 
s!- § = [ Yu do. (26) 


oO 


Squaring both sides of (26) and bearing in mind that 


= 2(1 - g's) = 2(1 - cos eVweeo 


(s' - s) 
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we find 


Ld Lo 
2 
= : 2 
€ f f V, V5 (44 He) dd, do, (27) 
Oo Oo 


Averaging the correlation coefficient for the refractive index fluctuations over the ensemble 


of realizations of the medium, we obtain 


Stet ey Lo LS 
C o 
e =u i [ VY5 N(x, - X52 ¥y- Yor 27 Zr) do, do, - (28) 
oO 


Since the curvature of the ray is small, the integration along the ray can be replaced by 


integration along a straight line: 
an ic, “Go Co 
eX = \a V_V. N(x.- x.) Ya- Yn» 2,- 2) dridr, . (29) 
12 J. ee ST, 2 2 Lo 2 
fo) fe) 
Introducing relative coordinates 


X= X- Xo VV] Vos 2 = 2,7 Zo) Y= T- Ly 


and center-of-mass coordinates 


se vee aa 
Re DP) I aa Volt 2 ee) ee Oe ee) 
we have 
V, VY, N(x, - Xos V7 Vos 247 Z,) = - TN ee) ; 


Since 45 >> a, we can integrate with respect to r from the limits - » to +. Eq. (29) takes 


the form 


= —_ +00 
Saat te f dr f 7 hea (40) 
O 00 
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Since the correlation coefficient N is an even function, we finally obtain 


—_ 


[oe] 
ea a as | S nar. 
J 
O 


If the medium is statistically isotropic and if consequently N depends only on 
"i 2 Ze: ee . ia : 
ee x +y+z , the Laplacian Should be expressed in spherical coordinates. 


then takes the form 


or 


where 


p=-30 [ gb, Os ee ON)" 7c. 
a 2 Or or 
O 


(31) 


Eq. (31) 


(32) 


55) 


(34) 


is a constant which characterizes the medium and is defined in terms of the correlation 


coefficient N. The quantity D plays the role of a ray diffusion coefficient. 


a 2 
Setting N=e 


a a 
Carrying out the integration in (34), we obtain 


2 


9 


Setting re =5X107, a = 60 cm, we find 


10 —-l 


. = ral 
D~1.5 X 10 af a Fc Se? an 


ai 


(35) 


-r/a ; 
If we use the other form of the correlation coefficient N =e / » then the integrand becomes 


i. © 2 ON 7 1 i-v/a pe. a7r/a 
2 or \* O&O - ae ar ; 


Because of the presence of the factor 1/r in the second term, the integral (34) will diverge 

logarithmically at zero. There is nothing unexpected about this. The correlation coefficient 
e- |F|/a corresponds to discontinuous changes in the refractive index fluctuation y, in which 

case the original ray differential equation (24) should not be used, since the gradient of 


refractive index fluctuations Vy becomes infinite at the points of discontinuity. 


5 The Angular Distribution of Rays. The Einstein-Fokker-Kolmogorov Equation. Since the 


properties of a medium with random inhomogeneities are given statistically, it is only pos- 
sible to predict the probability W(0,%,¢) that a ray which has traversed a path o will have 
the direction defined by the angles @ and g, if we are given the probability W (0,9; 0) for 
the initial point of the path (@ is the polar angle and d the azimuth). In other words, the 
propagation of a ray in a medium with random inhomogeneities can be regarded as a stochastic 
process without after-effect (continuous Markov chain), in which the role of time is played 
by the path o traversed by the ray. The distribution function W(9,0,6) satisfies the 


Kinstein-Fokker-Kolmogorov equation [25] 


ac. 3 2 or 
Qs Sea es QW) - 3g Ag QW) ee S66 QW) , (36) 
1 3 a 
+5 362 (Bug QW) + 30008 (Bog QW ) 
where 


g-sine, a- Bia B (37) 


oe (38) 
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The coefficients A and B are easily expressed in terms of the diffusion coefficient D, if we 


use the following formulas of spherical trigonometry: 


cos (@ + 48) = cos © cos € + sin @ sin € cos y, (39) 


sin Ag _ Sin yf 
sine sin (@ + Ae) ° (40) 


In these formulas ¥ is the angle of deviation of the ray from its initial direction. Since 


all directions of the deviations are equally probable 


sin f¥ = cos ff = O, sin“y = cos“ = 


KO) He 


Replacing the sines of small angles in (40) by the angles themselves, we obtain 


LO = sing 4(*in vy. (47) 
Moreover we have = — 
_— — o o 
—— € et e..|2OUCE ee ee 
= = 5 Sin y = 0, S37 Se (42) 
sin 9 sin @ 


Expanding (39) in powers of the small angles 4 and ¢€ and retaining in the expansion only 


terms of the second order in & and ¢€, we obtain 
2 2 
£6 +2 tan @ Ad + 2€ tan @ cos ¥- € =0. 
The solution of this equation has the form 


D 
48 = - tan @ + Vtan“@ = (2€ tan @ cos W- € ) 


(The second solution of the quadratic equation is not suitable since it gives a value for AO 
which is different from zero for € = 0.) Expanding the square root in powers of € up to the 


second power, we obtain 


A6 = - € cos ¥ +s ¢@ cot @ (1 - ae (43) 
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From this we get 


Moreover, it follows from (41) and (43) that 


Ap AB =0 . (45) 


Using Eqs. (33), (37), (38), (42), (44) and (45), we finally obtain the following formulas 


for the coefficients A and B: 


) 3 
(46) 
2D 
B.=2D, B ye Been 
ae op sine op 
Then Eq. (36) takes the form 
sin 9 Sy 2 (ein oH) yD TW : (47) 


which coincides with the equation for the rotational motion of Browmian particles [26]. For 


a ray starting out in the direction of the polar axis, Eq. (47) becomes 
; OW re) OW 
sin @ 55 = D35 (sin @ =) (48) 


The non-negative solution of Eq. (48), which satisfies the condition of normalization of 


probability and which reduces to zero everywhere for o = 0, 6 # O, has the form [27] 


W(6.@) = = » (2n + 1) P (cos 9) g7a(n41)Do , 


n=O 


where P_ (cos 9) are the Legendre polynomials. 
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We can form a more intuitive picture of the ray propagation if we calculate, for example, 


the mean value of the cosine of the angle of deviation of the ray from its initial direction: 


rl 

v= cos 9 = 2x . cos @ sin 9 Wd@ 
ws 
fe) 


Multiplying both sides of Eq. (48) by 2x cos © and integrating with respect to © from 0 to a, 
we obtain av/ao on the left hand side. The right hand side can be transformed in the follow- 


ing way by integrating by parts 


I se 
ex D [ cos 9 & (sin 9 ow } d@ = 2xD | cos 98 sin 9 a 
O O 
sig ne 
+ xD Fi sao a d@ = 2xD | sin 6iW 
e) O 
n1¢ 
- UD [ sin 8 cos OW dQ = - 2D . 
O 
For v we get the equation 
dv 
es 
whence we find 
v= soa s (49) 


At the initial point v=1. As ©—>o , v ~O, which has a simple physical meaning: after 
the ray has traversed a sufficiently long path 6, all of its directions become equally pro- 


bable. Similarly for the mean square of the cosine we obtain 


-6D6 
cos 8 = 5 (1 + 20 Be 


ol 


In the case of small deviation angles (sin 9~ 0) Eq. (48) takes the form 


Ow 1 W a 
i re 


The non-negative solution of this equation which satisfies the normalization condition has the 


Simple form 


al _-6°/lipe 


ae Lt DO 


i.e., for small deviation angles the angular distribution is Gaussian. For the mean square 


deviation angle we obtain 


of aie: . (51) 


Eq. (51) determines the fluctuation of the angle of arrival of the ray. Evaluating this 


fluctuation for a distance o = 1 km, and using Liebermann's data ( ue =5*x 10°”, a = 60 en, 


D~1.5 x 10° mt), we obtain \/ Of ~ 26.57 (minutes of arc). 


Kolchinski [28] analyzed data gathered from astronomical observations of fluctuations 
of the angle of arrival of rays from stars as a function of zenith distance. His analysis 


confirmed the relation (51) derived above. 


6. The Mean Ray Displacement. We now determine the mean Square rectilinear distance r 
from the initial point of a ray to the point at which it arrives after traversing a complicated 


path Oo in the medium. Using the formulas 


O 


O O 
x = ieee @d6o,ye [sin @ cos g dé ,25 [sin @ sin @ do, (52) 
J 

O O O 


2c 


we find 


ry 
II 
ra 
+ 
ca 
+ 
N 


(55) 
ome} 
= {| { cos 9, cos 9) + sin 8, Sin 9, cos a7 Py do, do,, 
0 Oo 


The expression in curly brackets is equal to cos €, where ¢« is the angle between the two 


directions (@,,0,) and (0,,9,). Therefore (53) takes the form 


The value of cos € can be obtained from (49). If Sn > o_, the deviation € from the initial 


direction (0, ,0,) is given by the formula 


-2D(6,-6, ) 
cos €,, =e 7 3 (55) 


If O5 < o> the direction (0,505) will be the initial direction, and the deviation from it 


is given by the formula 


cos € = e . (56) 


2) 


Taking this into consideration, we transform (54) as follows: 


on 
fo] C O al 
a) ees i 
r= f ao, [ cos €,, do, + f do, f cos Elo dG,, 
O Oo, O fe) 
Oo 6 an ) 
O -2D(6,- 3, ) l -2D( 0, - On ao 
2 do e 
= f do, f e do, + [ l 2 
O 
O Oo, O 
. ; 2D(6-6, ) a ; 2D6, 
- fmo-e Ja, + f w-e ) do, 
O O 
Carrying out the integration, we finally obtain 
2 6 1 -2D0 
Boe as ee ) (57) 
2D 
For the case of small Do (Do << 1) (57) becomes 
‘aoe 6 Ci = SN. 5 (58) 


which coincides with the formula obtained by Smoluchowski [29] for the mean square displace- 
ment of a heavy molecule which has traversed a complicated path o in a light gas. 
If the x exis is taken along the initial direction of the ray, the mean square ray dis- 


placement in the direction of this axis is given by the formula 
ome} 


a) Ee a eS 
x = fT] cos 8, cos 8, do, do, , (59) 


0 0 
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as follows from (52). We transform the right hand side of (59) as follows: 


6 o 6 of 
5) ie ee ee Re ee eae 
x ~f ao, f cos 0, cos e, d6,, +f do, cos @, cos eo, do, - (60) 
O 0, O O 


Eliminating cos @, in the first term by using the (already encountered) formula 


cos 8, = cos 9, cos E57 + sin eo, sin Eo, cos W 


P 


eliminating cos 9, in the second term by using the formula 


cos 9, = cos 9, cos E10 + sin 9, sin E15 cos ¥ >» 


and bearing in mind that cos vf = 0, we obtain 


O O Oo Oo, 
a) 2 
=f do, [/ c0s*e, cos €p cos “oe, cos E5 d6,, +f do, f cos 8, cos E15 do, . (61) 
O 7 O O 


; 2 
Designating the first double integral by Ij; we substitute for cos 8) and cos E57 using 


(50) and (55). This gives 


Oo Oo 
-6Do -2D(o.- o,) 
1 1 a =; 
rT, - | $0 + 2e )do, [ e do, (62) 
O O, 


2) 


Similarly, we evaluate I,, the second double integral in (61): 


. an -6Do. -2D(d_- o,) 
L, = do = (1 + 2e =) e ES Ge 
a ue 4 2 
O O 


Finally we obtain 


Eq. (64) takes the form 
‘aie & (1 - 2Do) 


in the case of small Do. We can now evaluate the mean square displacement of the ray 


initial direction. Using (57) and (64) we obtain 


‘ef = 
im Lea - e729) eS ie 
2D 18D 


In the case of small Do 


fs ae, 
V3 
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Using the "three-halves law", we can easily calculate the mean square displacement of a ray 


which has traversed a path 6 in the medium. Setting D = 1.5 x 10° jen we find 
V o% ~ 4.5m 


for a distance of 1 kn. 

Kharanen [5] also used the Binstein-Fokker-Kolmogorov equation in studying the problem 
of ray propagation in a medium with random inhomogeneities. He determined the ray deviation 
(angular and linear) from the plane containing the initial ray direction. It is easy to show 
that the differential equation obtained by him is valid only for small deviations and that it 
can be obtained from our differential equation in this limiting case. 

If we orient a spherical coordinate system so that the initial ray direction lies in the 
equatorial plane, then the ray direction with respect to the equator at any subsequent instant 
of time will be defined by the latitude X, connected with the polar angle 9 by the relation 
X+0Q=n/2. To find the distribution of the rays in latitude, we must integrate the dis- 


tribution function W(0,¢) with respect to azimuth @. Multiplying (47) by d@ and integrating, 


we obtain 


where 


Because of the periodicity of the function W the last term vanishes, and the equation takes 


the form 


t Pa t 
sino & = D cos 0 a eee ee 
9 
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If the deviation of the ray from the equatorial plane (9 = x /2) is small, then as an approxi- 


mation we can set sin 9 = 1 and cos 9 = 0. Then we get the equation 


Ow! | ow 
00 2 
OX 


which Kharanen used in his work. 


{. Fluctuations of the Transit Time and Intensity of the Rays. As conditions in the 
medium change, the time taken for a ray to go a given distance changes; moreover, the ray 
tube is deformed, which leads to intensity fluctuations. Let us consider the problem of 
transit time fluctuations and intensity fluctuations, assuming that the deviation of the rays 
from their initial direction (taken along the x axis) is small. 


According to (16) the time taken to go the distance L is given by the formula 


co 


Ole 


L 
f n(x,y, z)ax p) 
O 


O 


where the values of the refractive index n(x,y,z) are taken along the ray. The mean transit 


time is 


L L L 
Ree - 1 if it 
ee ee n(x,y,z)dx pee dx = u(x,y,z)dx, 
O O O 
O O O 


and finally the mean square transit time fluctuation is 


If L > a, then introducing relative coordinates x = Xo" Xo Y= Yo" Yyp 25 Zi Zz. we are 


justified in integrating with respect to x from the limits - © to» . Since N(x,y,z) is an 


even function, we obtain: 


nan oO 

ers 2 

ME =i I f N(x,y,z) dx 
“9 O 


The correlation coefficient N(x,y,z) differs from zero in the region where x has values of 
order a (x~a). Since the inclination of the ray is small, we have y << a and z << a in the 
region where the values of the correlation coefficient are appreciable. Therefore, in the 


last formula we can simply set y =z=Q. Then we obtain 


— “2 

i cut f N(x,0,0) dx. (68) 
Cc 
O O 


From this we can immediately find the phase fluctuation S‘' = wht: 


ge = Bie KL ; N(x,0,0) dx , (69) 
O 


where k is the wave number. 


ay 


The change of intensity is determined by the change of cross Section of the ray tube. The 


relative change of cross section as/S of the ray tube along the path dx is given by the formula 
as/S = div s dx. (70) 


Since the deviation of the ray from the x axis is small, we have 


3 Os Os. 
divs= ye + 3 ; Cia) 


and the ray equations (21) can be rewritten as follows: 


SS . Sx. » 
ax oy 2 
a(ns, ) he ; 
dx oe 


Integrating these equations, we find 


x 
x 2 Ob oy 
Lae > Sas 
O 
x 
zi Oh at 
“2 Oh [ oe. 
O 


Neglecting terms of the second order with respect to yw, we obtain 
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x 
t 
5 (4,952) -[ See axt 
O 


(72) 


xX 
Pe) ' 
s (x,y,z) -/ Sue ays2) axt 
2) 


Using (72) to replace s and s, in (70), we find 


era f Tu dx! , (73) 


where a designates the transverse Laplacian 


5° 


+ —— 


Pp. 
By? 5,2 


The relative change of intensity is connected with the relative change of cross section 


of the ray tube by the simple relation 


dL ds 
Tits: (74) 


Using (73) we find that the relative change of intensity along the path dx is 


; 2 


x 
SS ax [ at dx! 
@) 


To determine the change of intensity along the path L, we must integrate the last equation with 


respect to x between the limits O and L. Then we obtain 


pas 


L x 
log = = a ax | Fu ax! 
: O O 


x 
Setting u -/ wat dx' and dv = dx, we integrate (75) by parts. This gives 
O 


L 
log a4 7 | (L - x) wa dx. 


‘e) 


Replacing the intensity ratio by the ratio of squared amplitudes, we obtain 


NMlH 


log > = - 
O 


L 
f (L - x) Fi dx . 


We can use the last formula to determine the mean square amplitude fluctuations. 


(77) and averaging, we get 


x vi V5 Nyol%- Xr Vor Vys Zor 2) day dx 


Squaring 


Going over to relative coordinates x = X.~ Xs y= Var as: 2S Z~ zy and center of mass 


| a uit 4 | 
coordinates x = 5G X,) 5 Y= a(y,+ Yo) z= 5\Z,+ Ze) we can write 


x v v Nj o(x,¥,2) dx 
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in the case L>>a. Integrating with respect to x_, we obtain 
fe) 


i 00 
ae Ai N pol %¥s2) : (80) 


ey 4 / 


In the region of appreciable values of the correlation function we have X—~a, so that the 
2 2 
second term in brackets Lx /4t is small compared with the first term 17/3. Neglecting the 


second term, we get 


foe f)" - we i? f tv N ax 


Because of the smaliness of the inclination of the ray the integrand in this formula can be 


evaluated for y = z= 0. Thus we finally obtain 


(oe #) . : pe f ev? "| no (81) 


Eqs. (69) and (81) for the mean square phase fluctuation and the mean square amplitude 
fluctuation are not restricted by the requirement that the phase and amplitude fluctuations 
be small. In this respect they differ from the corresponding formulas obtained by Krasilnikov 
Gy who applied the method of small perturbations to the eikonal equation and to the equation 
expressing the conservation of energy flow in a ray tube. Using the same method of small 
perturbations, Bergmann [14 | studied the problem of phase and amplitude fluctuations in a 
spherical wave (point source). The formula obtained by him differs from (81) in that it has 
the full Laplacian under the integral instead of the transverse Laplacian, and the numerical 
factor is 1/15 instead of 1/6. 


Setting N = exp Rew os 22) /a2|, we obtain 


a 


22 


and 
(ee) 


gg le 


Eq. (81) takes the form 


2 —x 
8 a: 
(os £) Sees 
10) 


a 


If we set N = et/e, then 


ex] =a (fee) ov, 
y=z=0 
and the integral 
(ee) 
1 a -x/a 
[Ges ie 
5 ax ax 


diverges at the point x = 0. As already noted above, the correlation coefficient e7t/a 
corresponds to discontinuous changes of the refractive index, a case in which we cannot use 


the ray differential equation. 
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Part II 


DIFFRACTION THEORY 


Chapter III 


THE WAVE EQUATION 


8. Derivation of the Wave Equation for an Inhomogeneous Medium. The wave equation for 
an inhomogeneous medium differs from the wave equation for a homogeneous medium by the pre- 
sence of an additional term [55]. Accordingly, we do not think it amiss to give the deriva- 
tion of the wave equation for an inhomogeneous medium. We shall assume that the inhomogeneous 
liquid or gaseous medium is in a state of equilibrium. In the absence of the force of gravity 
(we neglect its effect) the condition for equilibrium is that the pressure be the same at all 
points of the medium. Thus, inhomogeneity of the medium will be caused by spatial change of 
the temperature and density of the medium; these changes cannot be regarded as independent, 
for, since the pressure is constant, a change in density is completely determined by a change 
in temperature. Strictly speaking, a medium which is inhomogeneous with respect to tempera- 
ture cannot be in a state of equilibrium, since temperature inhomogeneities are levelled out 
as a result of heat conduction. However, in weakly conducting media, the levelling out pro- 
ceeds slowly and heat conduction can be neglected in considering relatively rapid acoustic 
processes. Under actual conditions, in the ocean and atmosphere, the temperature inhomo- 
geneities are carried along by the flow and undergo convective motion due to the action of 
gravity. However, in the cases which interest us, the velocities of flow and convection are 
slow compared to the velocity of sound, and therefore the spatial distribution of inhomoge- 
neities can be regarded as quasi-static. The quasi-static condition will be discussed in 
detail in Section 26. 

Let Py denote the pressure, which is constant at all points of space, and let 
QO = p (x,y, 2) denote the density, which varies in space. Under the influence of an acoustic 


O 


wave, these quantities suffer changes and take the values p and p, respectively, so that 
— bd = 1) 
P i + Pl > P Pe + Pa , ( 


oy) 


where Py is the acoustic pressure and Py is the density change caused by the wave. The 


total pressure p, the total density p and the velocity of the acoustic oscillations v satisfy 


the hydrodynamical equations 


o & =v ; (2) 
oe = - div(pv) - (3) 


Substituting for p and op in these equations using (1), and regarding (as usual) the acoustic 
perturbations Py» Py and v as being small quantities of the first order, we discard terms 

of the second and higher orders of smallness. Then Eqs. (2) and (3) can be rewritten as 
follows: 


ee 
eS SR (4) 


Op) eee 
OL = "= iv(o v) . (5) 


Differentiating (5) with respect to time, we obtain 


apy are ov 
542 V\P 4 Ot : 


> 


since 0. does not depend on time. Substituting for O, a in the last equation using (4), 


we obtain 


3°, - 
2 = Vidas 2 (6) 
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To find one more relation between Py and Py» we must make an assumption about the thermo- 


dynamic character of the acoustic process. Considering the propagation of sound to be an 


adiabatic process, we can write 


3 


ap _ 2 do 
at at 


(7) 


where c is the adiabatic sound velocity at the given point of the medium. The equation con- 


tains total rather than partial time derivatives since the relation (7) must be satisfied 


for a given particle of the medium and not for a given point of space. Using (1) to substi- 


tute for the pressure p and the density p in (7), we obtain 


Si 
2 ot ~ ot“ yy ? 


Cc 


(8). 


with accuracy up to small quantities of the first order. In the case of a homogeneous medium, 


the convective term v-Yo. in the last equation vanishes, and therefore the relation (7) can 


be written for Py and o 


1 using either total or partial derivatives with respect to time. In 


an inhomogeneous medium the relation (7) can only be written using total derivatives. In going 


over to partial derivatives an additional convective term v-Vo. of the first order of smallness 


appears. The appearance of an extra term in the wave equation is the result of this convec- 


tive term. 


Differentiating (8) again with respect to time, we can write it in the form 


1 a 
— VS. Sg (9) 
S22 ot 0.62 2 


If we use (6) to substitute for the first term in the left hand side of (9) and (4) to 
substitute for dv/dt, we obtain 


3“ p 
di a ee : 
=— - Vp, + Vlogo,-Yp, = 0 . 


(10) 
a a+" 


3 


Thus, the additional term 


V loge. -VP, 


has appeared in the wave equation. Later we shall evaluate the order of magnitude of the 


additional term and explain the conditions under which it can be neglected. 


9. The Wave Equation for a Multi-Component Medium. In the preceding section we assumed 


that the medium is homogeneous in its composition. However, in practice, we can encounter 
the case where the composition of the medium changes from point to point, for example, ocean 
Water with variable salinity or air with variable water vapor concentration. In the equi- 
librium state, where the pressure is the same at all points, a spatial change of density is 
caused by a change of temperature and composition. In other words, the density is a function 


of the temperature and the concentrations of the separate components which make up the medium. 


We restrict ourselves to the simplest case of a medium consisting of two components. 
The generalization to the case of a multi-component medium is straightforward. Let o be the 
density of the medium, o' the density of the solvent, p'' the density of the dissolved 
component, and C the concentration, i.e. 
tf 
C=P>, p=o'sptt=pt(l +c). (1a) 
The density of the solvent satisfies the equation of continuity [30] 


Sp" = -div(p'v) , (12) 


where v is the velocity of the acoustic oscillations. A change in density of the dissolved 


component can be caused not only by convection but also by diffusion, i.e. 


Oo'! ap. eS 
R— = -div(o''v +t) , (13) 
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where 


F = ' = t 
I =-+ D,VC fe DVT : 


Here DS is the diffusion coefficient and D,, 1s the coefficient of thermal diffusion. Adding 


(12) and (13), we obtain 


3 i 5 
= = -div(pv) - divi , (14) 


for the total density o of the medium. In the absence of sound, this equation takes the form 


O x. ae 
as = -div I, ’ 


where I. is the flow caused by diffusion and thermal diffusion in the unperturbed medium, 
We shall neglect slow processes (and also the conduction, flow and convection of heat) and 
consider the distribution of inhomogeneities to be quasi-static, i.e. P4 to be independent 


> 


of time. Setting o = 0 (x,y,z) + ea =I. +i, in (14) and restricting ourselves to small 


quantities of the first order in Bay and T,, we obtain 


si = -div(o ¥) - div it, . (15) 


Finally, neglecting irreversible acoustic processes caused by diffusion and thermal 


diffusion, we write Eq. (15) in the form 


3p , 
—s = -div(o,¥) . (16) 


Moreover, neglecting viscosity and heat conductivity, respectively, we can write the two 


additional equations 


Ov 
Pe Ot = “VPy J (17) 


7 


dp _ .- a , (18) 


Thus, we obtain Eqs. (16), (17) and (18) for a multi-component medium; these equations 
coincide with Eqs. (4), (5) and (7) for a medium with a homogeneous composition. The only 
difference is that the parameters p. and a appearing in Eqs. (16), (17) and (18) depend not 
only on the temperature but also on the concentration at the given point. Repeating the 


considerations of Section 8, we obtain the wave equation (10) for the acoustic pressure. 
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Chapter IV 


SCATTERING BY INHOMOGENEITIES 


The theoretical study of wave propagation in an inhomogeneous medium reduces to the 
integration of a wave equation with variable coefficients, a problem of great mathematical 
difficulty that can be solved only for a few special cases. Here we must mention the exten- 


Sive investigations of Brekhovskikh [15255] on wave propagation in layered media, which he 


carried out over the last ten or fifteen years. The problem becomes simpler in the case of 

a weakly inhomogeneous medium, where the parameters appearing in the wave equation deviate 
only slightly from their mean values. In this case, the method of small perturbations can be 
successfully used. In connection with the problem of wave propagation in a medium with small 
inhomogeneities, the method of small perturbations has been used in two different forms: the 
usual form and the modified form of Rytov. Thus, for example, the usual form of the method 

of small perturbations was used by Ellison Ea Mintzer [21, 22,23], Tartarski [32] and others, 
in their work, while Rytov's method was used by Obukhov [12] and the author (47, 48]. Since 
in a certain sense the different formulations of the method of small perturbations complement 
each other, it seems to us expedient to study them both below and to evaluate their relative 
merits. For example, it is convenient to solve the scattering problem by the usual method, 
while it appears that it is more expedient to study amplitude and phase fluctuations in the 
direct wave using Rytov's method. We shall begin this chapter, which is devoted to scattering 


by inhomogeneities, by examining the usual form of the method of small perturbations. 


10. The Method of Small Perturbations. Dropping indices, we rewrite the wave equation 


(10) for the acoustic pressure 


2 
a op - va + Vlog e-Vp = 0. (19) 
as ato 


¥& 7 . 
Recently, further work (56,57, 61] using Rytov's method has appeared. In [58 | the classical 


method of small perturbations is used. 


var 


We assume that the density and sound velocity deviate only slightly from their mean values - 


and c_, i.e. 
O 


where 4y5 < oe and Ac < Coe Confining ourselves to quantities of the first order of smallness 


with respect to 4o and Ac, we rewrite Eq. (19) in the following form: 


2 2 
1 OD - OAc a p 1 
e© ate e a 8c 


Taking as the zeroth approximation the plane wave 


Py = A, exp [-i(wt - kx)] (22) 


(k = a/c. is the wave number in the medium with the averaged characteristics), we obtain the 


equation 


1 3p ee Ac op 1 
Se pe eo ee 5 - = Wp). Vp_, 
o° 52 1 2° 542 0 ( ) Pe (23) 


for the first approximation Py» or 


nae 


es Vp, = -|ax? == rae aad A, exp [-i(wt - kx)]. (2h) 
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Finally, introducing the notation 


_ [o2 2¢ , ik 3 | 
inQ = ax et 5. 2fo) A, exp[-i(at - kx} (25) 


O 


ho 


aS an abbreviation, we can write Eq. (24) in the form 

i 
5 > ee = 4nQ. (26) 
fe) 


Under the influence of the primary wave Ds each element of the inhomogeneous medium becomes a 
source of secondary scattered waves Py» where Q is the density of elementary sources. The 
total effect of the waves scattered by the volume V is given by the solution of the inhomo- 


geneous equation (26) 


Pp, = | ey , (27) 


where r is the distance from the scattering element (t,7,¢6) to the observation point (x,y,z), 


Les 


r= V(x-£)© + (y-n)® + (2-t)* 


Substituting the expression (25) for Q into (27), we obtain 


A 
2A ik 0 
note BEE 2A] 2 aot & o 
V 


(the factor e 1% is aropped). The scattering by fluctuations of the sound velocity is given 
by the first term in the square brackets, and the scattering by density fluctuations is given 


by the second term. 


11. The Scattering Formula. The fluctuations of velocity Ac and density 4&o which appear 


in Eq. (28) are random functions of the coordinates. The transition to ordinary functions is 
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achieved by squaring the expression (28) and then taking the statistical average. Assuming 
that the fluctuations of velocity and density are caused by fluctuations of temperature (for a 


constant pressure p,) we rewrite (28) as follows: 


A 
3 ° 
noo f PPS Gi, oo Ge, MEE] eolinten 


We shall calculate the pressure produced by the waves scattered by a cubical volume of side L 
which is considerably larger than the correlation distance a(L >> a) at a distance r which is 
considerably larger than the dimensions of the cube (r >> L). As can be seen from the diagram 


(Fig. 3) 


Fig. 3 The scattering configuration (schematic). 


ny is the scattering direction, n is the direction of the incident wave (coincident with the 
x axis), r' is the radius vector of the scattering point, r is the distance from the scatter- 


ing point to the observation point, and R is the distance from the origin of coordinates to 


the observation point. Then we have 


revR- nyr', (30) 


by 


obtain 


A | 
Py = = Ie | |x = ($2), ma a (So) - exp [ikR+ik(§-n, -r' av. (31) 
O O O O 


Before proceeding, we note that Eqs. (30) and (31) can be used only in the Fraunhofer zone. 
In fact, Eq. (30) is obtained as a result of expanding the exact expression for r in a series 
of powers of r'/R and dropping all terms of the series beginning with the quadratic terms. 
This means that the contribution to the phase produced by the first discarded term (the 


quadratic term) must be small compared to nm. Expanding the expression 


Since the largest value of r'~L, the condition that the quadratic terms be negligible can be 


written in the following way: 
2 
kL" /R << x. 


This condition is satisfied for the Fraunhofer zone. The minimum dimensions of the scattering 
volume are of the order of the correlation distance (L~a). Consequently, we obtain the 


least restrictive condition for R by setting L = a in the preceding inequality. This gives 


ca" /R ale 


Multiplying (31) by its complex conjugate, we obtain 


A5 


A 4 ceeS 5 
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fe) 
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x exp |i (eres) - ee - =] dv, dv, . 


«| (84 -B) Ba Grt 2) | = fay - Th) - nyt} - 2) (33) 


Where we have introduced the notation 


where @ is the scattering angle; r is the distance between points inside the scattering volume. 
Taking the statistical average of (32) and introducing the correlation coefficient for the 


temperature fluctuations 
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(an)* 
we obtain 
2 O 
A’ (AT) 4 D) 
2 “0 4k" ,0ce\2 2ik- ,d yo 3 3 
pl -2o5 ff MS ge? we SES Le Sa 
1 ESE oT (so)... ) (34) 
16x R vy ‘Y ou : Py Po t Po %o or Py 38) 3e,, 
2 C >> 
Kk sOpne ON iKer 
+= (+) Ne oe e dv, dv 
52 Dp, £1085 1 e 
Now introducing the relative coordinates 
f= 6,7 §5 
N= 1, 7 No» (35) 
¢ C1 7 bo» 
and the center of mass coordinates 
1 
J. 
ule 
z= 5 (6, + $,), 


and bearing in mind that the correlation coefficient N(r) depends only on the relative 


coordinates, we calculate the derivatives of the correlation coefficient which appear in (54) 


SN | SN oN aN an | ON 
de, 8B 7 EE 96,06, ee 


uf 


Then Eq. (34) takes the form 


2 2 
AW (AT) 4 oes, 
2 O uk” ,0e,\2 Pike ,O 1 ,dc ON 
p= Soe ff AS Gps n+ SE Gh, & Gap 2 Se (37) 
V O 
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The integrand does not depend on the center of mass coordinates. Carrying out the integration 


with respect to these coordinates, we simply obtain the factor V. The integral with respect 


to the relative coordinates reduces to three integrals, and integration by parts is used to 


evaluate the integrals containing derivatives of the correlation coefficient N(r). In fact, 


=> 


noting that Ker = KLE + Kan + Kos Wwe obtain 


E 
aie or a dgdnde Sf Es A dndt - if iK, eth'T vaeanat 
1/2 -L,/ 2 
= 20K, f ee tay: 4 
V 
+L/2 0, +L/2 tL /2 +L/2 2 
[BF 2 was = PP EEF as FPP age 8 ee 
-1,/2 s iP aie -1,/2 


of the cube, 


mitt 


since the correlation coefficient and its derivative vanish on the faces & = + 


the dimensions of which are large compared to the correlation distance. Because of these 
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relations, Eq. (37) becomes 


2 2 2 
A” (AT) 2 K,k ae 
2 O 2k” 0c 3 > ; 
Ip, [P= v = Gip +t GB, | f mG) oF a 8 
where 
K, = Ken = k(A - BL) R= (1 - BL-R) = K(2 - cos 0) = 2k sin? 8 (59) 


Using (39) to replace K, in (38), we finally obtain the following formula: 


E 


A” (AT) & 2 ee 
Ip, |° = ag le (SS), + = ($2) sin® a f N(r) e> ” av. (40) 
x O 
1 


Eq. (40) lets us compare the scattering by velocity fluctuations with the scattering by density 
Fluctuations: these correspond to the first and second terms in the square brackets preceding: 
the integral, respectively. 

In the case of a gaSeous medium at the constant pressure Po the sound velocity fluctu- 


ations are given in terms of the density fluctuations, as can be seen from the formula 


a = ass 
le } 


for the velocity of sound, so that 


= (Sp 7 Gy (41) 


O O O 


Consequently, in a gas the scattering by velocity fluctuations and the scattering by density 
fluctuations agree in order of magnitude. In liquids Eq. (41) is no longer valid and the 


left hand side predominates. Thus, for example, according to the data of Beranek [a4], 
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for sea water at 15° C and with a salinity of 36%, i.e., the velocity fluctuations are an 
order of magnitude larger than the density fluctuations. For this reason, bearing in mind 
subsequent applications of the scattering formula in hydroacoustics, we shall neglect density 
fluctuations in comparison with sound velocity fluctuations. This cannot usually be done in 


atmospheric acoustics (this matter is discussed in more detail below). For water Eq. (40) 


simplifies to 


Gr. J N(F) ev’? av, ie) 
_ oe) 


Note that neglecting density fluctuations is equivalent to neglecting the additional term 
V log p-Vp in the original wave equation (19). 

The form of Eq. (42) can be abbreviated a bit by introducing the refractive index n in- 
stead of the velocity ec (n = c./e). Denoting by p the deviation of the refractive index from 


its mean value (equal to unity), we find from n = ie that 


PO et Ge Ney AR 
O 0 
>) 
= 2 
ue =45 ($8) (ar)* . 
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Eq. (42) can then be rewritten in the following way: 
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IP, | = V eer rf N(r) as v av. (44) 


This is the form in which it is encountered in Pekeris' paper [S5 ls 

If the statistical properties of the medium are isotropic, then the correlation coefficient 
depends only on the modulus of r. In this case it is expedient to introduce a spherical 
coordinate system, the polar axis of which coincides with the vector K, and to carry out the 


angular integrations, i.e. 


o 1 On 
ikr cosa 2 
e r 


N(r) dr sin adad¢ - i [ N(r) sin(Kr) rdr. 
O 


O00 


(The upper limit in the integral with respect to r can be set equal to ~, since the dimensions 


of the cube are large compared to the correlation distance.) Eq. (44) becomes 


ee a 00 
a A, k LU 
P| — an Sar ean f N(r) sin(Kr) rar. (45) 
mR K - 


The integral in this formula can be evaluated by specifying the form of the function N(r). 


For example, setting 


Nr) = eT, (46) 


we obtain the tabulated integral 


- 3 
f et/2 sin(Kr) rdr = — 5 
O ( + gar te : 
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In this case the scattering formula can be written in the following way: 


a. eB 
ic 
eal Os 7 (47) 
Pyl=¥ e Ba. eee 
mR (2+ mate sin 3) 


This formula can also be found in the paper of Pekeris cited above. Setting 
-r* /a® 
N(r) =e 3 (48) 


we again obtain a tabulated integral and the scattering formula becomes 


Ie 3, Be 28 

Se ee a6 2:2. BE 

ee exp|-k a sin 3 : (49) 
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From Eqs. (47) and (49) we can draw the following conclusions, which will be important later: 

1. "Small-scale" fluctuations (ka << 1) are the cause of isotropic scattering. 

2. In the case of "large-scale" fluctuations (ka >> 1), the scattering has a sharply 
directed character, and most of the scattered power is concentrated within the small solid 
angle © ~1/ka. 

In atmospheric acoustics it is usually necessary to take into account density fluctuations 
as well as velocity fluctuations. As can be seen from Eq. (40), the scattering by density 


fluctuations corresponds to the term 


which has a maximum in the direction 9 = xn (reflection). However, in the case of large-scale 


inhomogeneities this term is suppressed by the factor 


fue Pe 


Vv 
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Which has a strong minimum in the same direction. In other words, in atmospheric acoustics 
we can neglect density fluctuations as compared with velocity fluctuations, for large-scale 
inhomogeneities; in hydroacoustics we can do so for any relation between the wave length and 
the size of the inhomogeneities Moreover, these conclusions remain valid in the case where 
the velocity and density fluctuations are caused by concentration fluctuations. In this case 
the formulas will contain derivatives with respect to concentration instead of derivatives 


with respect to temperature. For sea water [30] We have 


ay Oc = ) 
a (So), p, ~0.8 , D (30)0,p_ ~ 0.15 - 0.20, 


i.e., the velocity fluctuations are 4 to 5 times larger than the density fluctuations. 

We can obtain similar scattering formulas for electromagnetic waves. They differ from the 
formulas obtained above using the scalar wave equation by the presence of the factor sin’X, 
Where X is the angle between the scattering direction and the direction of the electric vector 
in the incident wave. This factor takes into account the nature of the polarization of the 
incident Wave. 

In their theoretical investigation of the scattering of radio waves in the troposphere 
Booker and Gordon [36] assume that the correlation coefficient of the dielectric constant 
has the form et/2, as a result, they obtain the scattering formula (47) which is valid in 


oa 2 


this special case. Fejer [42 | used a correlation coefficient of the forme” and obtained 


Eq. (49). Staras cd applied the method of small perturbations to Maxwell's equations with- 
out any special assumptions about the correlation coefficient and obtained the general scatter- 


ing formula (44). 


12. The Scattering Coefficient. The attenuation of energy flow in a wave which passes 


through a layer of thickness L is given by the formula AI = alIL, from which we obtain the 


relation 


a - ot (50) 


for the scattering coefficient. The energy flow I through the boundary of the cube is pro- 
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I= pA TL, (51) 


where B is a constant of proportionality. The attenuation AI of the flow in going a distance 


L is equal to the total energy scattered by the cube, i.e. 


AI -p Ip, | as. (52) 


Here the integral is taken over the surface of a sphere of radius R. Confining ourselves to 


the case of a statistically isotropic medium and using Eqs. (45) and (50), we obtain 


a) ~ 
a = <i z dn [ N(x) sinlte) rdr, (53) 
O 


where dQ is the element of solid angle. Recalling that K = 2k sin © and changing the order 


2 
of integration, we rewrite Eq. (53) in the following form: 
ake oe) 14 
Yo 
a=2u 1K? f N(r) rdr Fi sin(2kr sin 3) cos = de. (54) 
O O 
Integrating with respect to 0, we finally get the formula 
2,2 
G S220 2 [1 - cos(2kr)}] N(r)dr. (55) 
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We now calculate the scattering coefficients for the different correlation coefficients. 


Setting N = er /2 we find 


54. 


_— 8 a= a? 
1+ 4k°a 
Pea. 
Setting N=e” /* | we obtain 
— oe. 
a= Sx Ve k-a(1 sere ). (57) 


In the case of small-scale inhomogeneities (ka << 1) these formulas become 
-@ dt 20 
a=8u° k a’, a= \/x i te, (58) 


Thus, as was to be expected, for small-scale fluctuations we get the Rayleigh dependence: the 
scattering coefficient is proportional to the fourth power of the frequency. For large-scale 


fluctuations (ka >> 1), we have 
Be 2 2 
a=2u ka, a= Jn um ka, (59) 


i.e., in the case of large-scale inhomogeneities the scattering coefficient increases as the 


square of the frequency, like the coefficient of viscous absorption. 


13. Applicability of the Scattering Formula. The condition for the applicability of the 


scattering formula can be obtained from the requirement that the scattered energy be small 


compared to the incident energy, i.e. from the condition AI << I, which by (50) gives 
aL << l. (61) 


Considering the case of large-scale inhomogeneities, we use (59) to replace a in (61) and 


obtain 


eee, (62) 
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In particular, if the scattering volume has dimensions of the order of the scale of the 
inhomogeneities (L ~a), then the condition (62) for the applicability of the scattering 


formula can be written as 
the Ka’ << 1. oor \/ ie ka << l. (63) 


The inequality (63) also follows from an analysis of the suitability of an approximate solution 
of the problem of scattering by a sphere of radius a, whose properties deviate only slightly 
from the properties of the surrounding medium Bae In our case, of course, the role of the 
radius is played by the characteristic scale of the inhomogeneities (the correlation distance). 
It follows from the condition (62) that the dimensions of the scattering volume and the 
scale of the inhomogeneities must not be too large, and the frequencies must not be too 
high. In particular, this means that it is impossible to make the transition to geometrical 
acoustics (k > o) in the scattering formula. In this case the formulas of the preceding 


Section would give infinite values. 


14. Attenuation of an Acoustic Bundle Due to Scattering. The results of Section le 


pertain to a plane-parallel bundle. If the bundle diverges, as is often the case, then in 
calculating the scattering coefficient it is necessary to take into account that part of the 
energy which is scattered out of the bundle. The energy scattered inside the bundle does not 
cause it to be attenuated but rather causes another effect, namely amplitude and phase 
fluctuations in the direct bundle (see Chapter 5). 

If an axially symmetric bundle has a divergence angle 28 os then in calculating the 
scattered energy and the scattering coefficient using Eq. (54) it is necessary to integrate 


With respect to 6 from the limits oe to x. Then we obtain 


ss i 8 
a= 2 fe i f n(x) [eos (2 sin 2) - cos(2kr | dr. (64) 
O 


a 
Setting N=e /2 and carrying out the integration in (64), we obtain 
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>, a sin“ > 
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If N=e » then (64) gives 
= 9 22 
a= \/x < rae vxp| x2e2stn2 = eee (66) 


In the case of large-scale inhomogeneities most of the scattered energy is concentrated inside 
the angle @ = 1/ka. This angle can be larger than or smaller than the divergence angle 0. 
of the bundle. If @< Q,»most of the energy is scattered inside the bundle. In the other 
limiting case (6 >> 8.) most of the scattered energy leaves the boundaries of the bundle. In 
the latter case in calculating the scattering coefficient we can use the formulas derived for 
a plane wave. 

Liebermann [4] gives an incorrect formula for a plane-parallel bundle, which differs from 
Eq. (56) by the factor eae as a result Liebermann gets a value of the scattering coefficient 


which is three orders of magnitude too small: In fact, using his experimental data is = 


5 x 10°”, a = 60 cm) Liebermann gets a ia dore en Pere frequency v = 10 kilocycles. The 


calculation based on Eq. (56) gives a = 10! ain = 9904 db/km, which is approximately a 
twentieth of the observed attenuation in sea water at the same frequency. Since the scatter- 
ing and absorption coefficients depend in the Same way on frequency (proportional to the 


square of the frequency) the ratio just obtained between the scattering and absorption 


coefficients remains the same (in order of magnitude) at other frequencies too. 


Chapter V 


FLUCTUATIONS 


When a wave propagates in a medium with random inhomogeneities, fluctuations of the 
characteristics of the wave field due to the superposition of the scattered waves and the 
primary wave are observed. There must be a dependence between the fluctuations of the 
characteristics of the wave field and the fluctuations of refractive index. The theoretical 
problem consists in finding this dependence, which can then be used to draw conclusions about 
the statistical properties of the wave field, given the statistical properties of the medium. 
The inverse problem is usually not unique and cannot be solved without additional assumptions. 
However, if we make reasonable assumptions about the form of the correlation coefficient of 
the refractive index, for example, then by measuring the field fluctuations we can determine 
the mean value of the refractive index fluctuations. Thus, the study of wave fluctuations, 
which is interesting in its own right, opens new possibilities for the study of the properties 
of the medium which the wave traverses. 

We begin our study of the statistical properties of the wave field by examining amplitude 
and phase fluctuations. In order to do so, we shall derive the basic formulas for the ampli- 
tude and phase fluctuations by two different methods; this will allow us to compare the merits 


of the two methods. 


15. The Method of Small Perturbations. In order to be explicit, we shall assume that the 


random inhomogeneities occur only in the right half-space (x > 0) and that the left half-space 
(x < 0) contains no random inhomogeneities. A plane wave 
16, -i (wt-kx) 


y= A. e = A. e : 


advances from the homogeneous to the inhomogeneous medium. A receiver is located in the 
inhomogeneous medium at the point with coordinates X,y,Z, and the waves scattered by the 


innomogeneities as well as the inci c . . 
wave Py are incident on the receiver. We shall describe their 
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total effect by the wave function 


ig, 
3 


ali 


which can easily be found to a first approximation by using the method of small perturbations, 
if we take D. as the zeroth approximation. Neglecting density fluctuations as compared with 


velocity fluctuations and using (28), we obtain 


k A, aik(r+é ) 
f ra aa BCE yase) dv , (67) 


where r is the distance from the scattering element dv with coordinates €,n,¢€ to the observa- 
tion point (x,y,z). The integration in (67) is to be taken over that part of space from which 
scattered waves arrive at the observation point. If we denote the result of superimposing the 


primary and scattered waves by 


Dividing both sides by Ae °, we obtain 


i($-6.) 44 1 ity Po) (68) 


>| 
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where by (67) the last term is given by the formula 


a - 2 ; 
2 e cee ES J 2 gtelr-(e)] u(E,7,6) av . (69) 
- V 
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Introducing the notation 


“ AP)! er (69a) 
we find 
i ab 
X = os f 7 cos k[r - (x - é) | Lesa) “ay, (70) 
ie dk 
te On [ - sin k[r ~ x = é) | Ley 955) dv . (70a) 


Since the amplitude of the scattered waves A, is small compared to the amplitude of the pri- 
mary wave Aw the amplitude A = A. + OA and phase gp = p+ xt) of the resulting wave differ only 
slightly from the amplitude A, and phase p. of the primary wave, i.e. 


— << 1, Ap <<. 


Expanding the left hand side of Eq. (68) in powers of the small fluctuations and keeping only 


leading terms, we obtain 
~ + iA@=X + iY. 
O 


Equating real and imaginary parts separately and using (70) and (70a), we obtain the following 


formulas 


2 
dp = < : = sin k[r - (x - &)] u(,n,6) av, Gan 
AA Ke a 
AS anes f r cos k[r - (x = E)] PCE snc) dv (72) 
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for the phase fluctuation A and the amplitude fluctuation AA. In these formulas the argument 
[r-(x-€ )] gives the difference in path length to the observation point (x,y,z) traversed by 
the primary wave and by the wave scattered by the element Cena Cs We note that this differ- 
ence becomes zero if the scattering element lies on a straight line through the observation 
point and parallel to the x-axis, as is to be expected in this case. 

The expressions just obtained for the amplitude and phase fluctuations are restricted by 
the requirement that I) << 1 and dA/A.. << 1, and consequently by the requirement that the dis- 
tances be small, since the fluctuations increase with distance. Therefore, it is more expedi- 
ent to use the method of small perturbations in the form given by Rytov [11] in connection 
with the theory of diffraction of light by an ultrasonic lattice. As we shall see below, in 
this form the method of small perturbations is not subject to such a Stringent restriction; 
this is essential for comparing theory with experiment, since experimentally one often 


observes large amplitude and phase fluctuations. 


16. Rytov's Method. The acoustic pressure satisfies the wave equation 


2 2 
as - LL) op a Te = 0, (74) 
fas st2 
Co 


where w is the deviation of the refractive index of the medium from its mean value of unity 


and |u| << 1. As before, we write the plane wave in the left half-space in the form 


in eh e7i(wt-kx) - 
O O 


We are looking for a solution of the wave equation in the right half-space in the form 


p = A(r) _i[wt-s(r )] ; (75) 


where A(r) and. S(r) are unknown functions. The gist of the method consists in replacing the 
wave function p by another function y which is connected with the first function by the 


relation 
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a et lot-W(7)] | (76) 


= fo) 


It follows by comparing (75) and (76) that the function y(r) is given by the following relation: 


V(r) = S(F) - i log A(z) ; (77) 


O 


As is apparent, the real and imaginary parts of the function y(r) just introduced determine 
the phase and the logarithm of the ratio of amplitudes (sound intensities), respectively; the 
latter are the fluctuations in which we are primarily interested. 

Using (76) to replace the function p by W in Eq. (74+), we obtain the following equation for 


the function yp: 

(Vy)> - ay = ken’. (78) 
The zeroth approximation ¥ satisfies the wave equation for the homogeneous medium, i.e. 

(vy.)* J ivy i. (79) 
Subtracting (79) from (78), we obtain the equation 

a(vy Vy") - ivy? = 2Quk> + es - (4) (80) 


for the function yw' = y - Vo Using the method of small perturbations and assuming that Vy' 
(more accurately, the dimensionless quantity = Vy") is of order yw, we drop the terms contained 
in the square brackets, since they are small quantities of the second order with respect to uh. 
Then for y' we obtain the linear equation 


a(Vy Vy!) - ay" = auKe, (81) 
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This equation was obtained under the assumption that = Vy > A pw or 


1 t 
k | Vy | << (82) 


This last condition means that the phase change and the relative change of amplitude per wave- 
length must be small. The relation (82) imposes no restrictions on the total change of these 


quantities. Setting ¥, = kx in (81), we obtain 
! 
: 2 
2k oe - ivy! = 2uK°. (33) 
Introducing the new function W defined by 


yi=e W, (84) 
we find for it the equation 


Siew = f(xy yz); (85) 


where 


f(x,y,z) = 12uKe a. (86) 


The solution of the inhomogeneous equation (85) has the form 


Wa- [ =e £(£,n,6) av, (87) 


where 


rey (x-t)© + (y-n)* + (z-t)® . 


We must take the plus sign in the exponent so that the formulas of this section will reduce 
in the case of small amplitude and phase fluctuations to the formulas of the preceding section, 


which were obtained by using the usual form of the method of small perturbations. 
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Using (86) to replace f, we obtain 


ee ; 
as i Da ae) can, (88) 


ag 


Using (84), we now return to the function y' 


gee f 1 ikle-(-8)} ce nt) av. (89) 
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Using (77) to separate real and imaginary parts, we obtain the following formulas 
oa lh 
s'=- => sin k[r - (x - é) | udv , (90) 


—> 2 
log AW) - KU i cos k[r - (x - é)] udv (91) 
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for the amplitude and phase fluctuations. In the case of small fluctuations they reduce to 
Eqs. (71) and (72) of the preceding section, respectively, which were obtained by using the 
usual form of the method of small perturbations. 

We now discuss the limits of applicability of Rytov's method. These limits are a conse- 


quence of the inequality (82) 
ile 
= [wi] «<1. (92) 


This inequality implies the smallness of the relative amplitude change over a wavelength 


<< (93) 
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and the smallness of the change of the phase fluctuation over a wavelength 
iL t 
. M'| < i. (9) 
The inequality (93) is met if the amount of scattering in going a wavelength is small. 


The inequality (94) implies that the angle of inclination of the ray to the initial direction 


is small. In fact, since 


S = 5. +S’ = kx + S'(x,y,z), (95) 
we have 

oS _ os! os as! oS os! 

ae =k + a ) Oy = Sy ) Ae = Aa . (96) 


By (94) the relations 


0s 


Selo 2 


37| <<k; (97) 


are valid, 1.e., the transverse components of the gradient of phase are small compared to the 
longitudinal component. Therefore, the angle of inclination of the ray with the initial 
direction is small. 

The condition that the amount of scattering in going a wavelength be small is always ful- 
filled in a weakly inhomogeneous medium (lu| << 1). Since large-scale inhomogeneities produce 
Sharply directed scattering, in this case the condition that the angle of inclination of the 
ray be small is also met. Small-scale inhomogeneities produce isotropic scattering, so that 
in this case the condition that the angle of inclination of the ray be small is met only when 


the amplitude of the scattered waves is small compared to the amplitude of the incident wave. 


Laie Comparison of the Methods. After studying both methods and finding their limits of 
applicability, it is natural to compare their relative merits. The usual form of the method 
of small perturbations is subject to a distance limitation. In the case of large-scale inhomo- 


geneities, Rytov's method is not subject to a distance limitation, and is therefore more gene- 


65 


ral. In the case of small-scale inhomogeneities the two methods are equivalent. Since large- 


scale inhomogeneities are observed under actual conditions, we shall give preference to Rytov's 


method. 


18. The Fresnel Approximation. In the case of large-scale inhomogeneities (ka >>1) we 
can neglect wave reflection and limit the region of integration in Eqs. (90) and (91) to the 
layer which lies in front of the receiver and is contained between the planes & = O and 
€ = x. An appreciable effect will be produced only by those inhomogeneities which are con- 
centrated within a cone with its vertex at the receiving point and with an aperture angle of 


the order of 1/ka. Inside this cone the formula 


r=V (x-E)° + ae ; 


where 


po SG aay se 2), 


can be replaced by the approximation 


1 C 
ra(x-e) +3 22). (98) 


Using (98) to substitute for r-(x-&) in (90) and (91) and replacing 1/r by 1/(x-&), we obtain 


a ne sin [kp /2(x-t J] 
cee Te 
S' = se fff en aad u(—E,n,6) d&dnde, (99) 
O = & 


. 2 2 COs [ico /2(x-8)] 
tog ADE [Tf ; u(E,ns6) dana. (100) 


A 
O x 
oO - & 


Eqs. (99) and (100) correspond to the Fresnel approximation in diffraction theory. In 
Appendix I we give a more rigorous justification of these equations, using the Sommerfeld 


integral. 
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We introduce the following abbreviations in the integrands of (99) and (100): 


oe 
x 
D 


6. = 6) e Sin ) 
1\ "kK? Ot Caan: B(x - €) 


(101) 
0.6 oe = ss cos io" ‘ 
2 k ? en(x - & Mecca £) 
Moreover, we drop the prime in writing the phase fluctuation S' and denote the amplitude 
fluctuation by B = log A/A.. Then Eqs. (99) and (100) can be rewritten as 
xX + 7 
xX ~ 
S(x,y, z) =k fff $, ( xe »P ) NCE ns) déande , (102) 
O - © 
X +o 
xXx = 
B(x,y,z) =k [ff $, ( xa , ©) wlé,n,$) dédndg. (103) 
O - @ 


Eqs. (102) and (103) were obtained by Obukhov [12]. We now introduce the dimensionless vari- 
ables 
x' = kx, y' = ky, z' =kz, €' = ke, n' = kn, €' = ko, p’ = kp. 


Then finally Eqs. (102) and (103) can be written in the following form: 


x'+ 0 


[LP Bot - ee) wean asranage, (10!) 


O - © 


UI 


SCX sy a2) 


x'+ 


PEF Boost ~ ste) wets nts8") agtantage (105) 


Oo - © 


lI 


BCX! jo 52") 
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19. Amplitude and Phase Fluctuations. Turning now to the statistical aspects of the pro- 
blem, we undertake the calculation of the mean square amplitude and phase fluctuations. 
We assume that the receiver is located at the point (L,0,0) or (L',0,0), where L' = kL. Then 


Eqs. (104) and (105) can be rewritten as 


L'+ © 
S(L',0,0) -f[f 6, (E'-€',0') wets n's6')ag'an'ag', (106) 
O- w 
L'+ © 
B(L',0,0) -[ ff $,(L'-§',0') Le sh ode tan dc’: (107) 
O - 


2 2 
where p' =Vn' + €'  . Squaring (106) and (107) and averaging, we obtain 


L'L' +o 
e) 2 t t t t t , 
S =u JILL LL §, (£'-£4, 61), (L -§5,05)N(r' dg Idésdnidnsacs acs, (108) 
OO = 
L'L' + 0 


NM 


Baw PLL LLL Beep eB co-sy, op w(x" passat santansacraty, (109) 
00 -o@ 


where N(r') is the correlation coefficient of the refractive index and 


(e1- 68)+ (nt= ng) (tr 63). 


We shall consider only statistically isotropic media. In this case the correlation coeffi- 


client N depends only on the modulus of r. Introducing relative coordinates 


Vem - No $= 6 - ob (110) 
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and center of mass coordinates 


y=Z(nt ens), 223 (tr + th), (111) 
we obtain 
—  _ LL +2 
Fe PPP [f a be -u Varn? + Gea? (212) 
OO —- © 


x >, ee CG ae (2 - 2)? | N(r')d&5 dé tdnd€dyaz. 


L'L' +o 


ee TTT TE bbe, Jaen? Gra? ] (113) 
x 9, [ ae (4 a (2 - 2 | N(x! )dgsd&sdndfdyaz. 


We can carry out the integration with respect to the variables y and z in these formulas. 


From the evaluation of integrals of similar form given in Appendix II, it follows that 


+ ow 


I , Gi 8p vg e y's G+ a | (114) 


- © 


x $, cE = Es VG z y)* + (2 : OF foyas 


2 Gtey - ey 0) Bla - ey + eps a) 
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+ 
ffab-w Gore): os 


where - = n° ee aor 


After the integration, Eqs. (112) and (113) become 


BF 5 te oe + Ty); (116) 
B = : ue (I, - I,), (117) 
where 
L'L!' + © 
er eat $,(é - §3, 0) M(r')ds jag sande, (118) 
0OO- @ 
LiL L* += 
Ty “If Re ee ar) Gay ole ck ler (119) 
00 - @ 


Further simplification of the integrals IL and I, 


is large compared to the correlation distance (L >>a). Introducing the relative coordinate 


is possible in the case where the distance 


a EF - BD and the center of mass coordinate x = a(t! + EA), we are justified in this case in 


integrating with respect to € from the limits - ~ to +» . Then we obtain 


+ 00 
wx f ff 9, (6,p) N(x! )dganas, (120) 


{O 


Ls 


sf 


O 


L! 00 
I, ay se [PP Bev - 2x, op) N(r')dédndc. (121) 


Since the integrand in (120) does not depend on x, we immediately get 


A ffl 9, (€,0) N(r')aganae (122) 


FH 
tf 


2 
a ; 
uv fff Bae sin or N(r')d&dndc. 


In (121) it is also possible to carry out the integration with respect to x. Since only one 
factor in the integrand, namely $, [eu - 2x, |, depends on x, the problem reduces to calcu- 


lating the integral 
L! 


L! ; 
iL. , QO 
t - a 
fi & 2x, 9] dx J an(OLi- ex) 1" BUBET- oxy 
O 


O 


By the substitution 


2 


Ae? “DER By) 


it is reduced to the form 


a 2 
al: sin z cx op, al : Q 
i [ sine ag == de at (£7 


o/b" 


(see the definition of the integral sine (39]). Thus 


if a 


Ee [ff a si (ence depanat. (123) 


If in (122) and (123) we go over to polar coordinates (p,@) in the (n,¢) plane, and if we bear 


in mind that the correlation coefficient N(r') is an even function of &, we obtain 


[o.6) [o.6) 
1 se 
= t _— i —_— ’ 
I, = eb [ as | f Sin 3 N(r') pdp, (124) 
O O 
co oo 
= Pp 
a --[ at | si (7 N(r' )pdp (125) 
O O 
Eq. (124) can be simplified further. Instead of po we introduce the new variable 
2 
q =p /2t. Then the second integral in (124) can be written as follows: 
co 
Fi sin q-N(r')dq. (126) 
O 
We now integrate (126) by parts twice, i.e. 
ioe) o6) [o.6) 
: S , ON(r') 
sin q:N(r')dq = | - cos q-N(r') + cos q zi dq (127) 


O O oO 


oo 


N(E,0,0) + 


i 


oO D 
sin q - -/ sin q SMe") dq 
oq 
O 


O 


[eo «7 > 
N(E,0,0) -/ sin Q 2 dq. 
oq 
O 


fe 


As p changes, the correlation coefficient falls off from its maximum value of N(E,0,0) at 

o = O to values which are effectively zero at distances p of the order of the correlation dis- 
tance a'(a'= ka). Correspondingly, the variable q changes to a quantity of the order of mag- 
nitude aoe. and since the region of appreciable values of & also does not exceed a‘, the 


change in q will not be less that a'/2, Lees 


In view of this we have 
: (129) 
ne 


: il: 
Therefore, the integral in the right hand side of (127) is of order —rN(E 0,0). In the case 
under consideration, namely large-scale inhomogeneities (a! >> 1), we can neglect this inte- 


gral as compared with the first term N(€,0,0). Then we obtain 


co 


[ sin q¢N(r')dq = N(E,0,0). (130) 


O 
The integral I, (122) reduces to the form 


oo 


i= aut | N(E,0,0)dE. C152) 


O 


By introducing the variable 


+ 26 132) 
ome Te (13 


the integral I, (123) reduces to the form 


ae aut f at | si veN(r')dv. (133) 
O O 


(3 


Finally, for the mean square amplitude and phase fluctuations we get the following formulas: 


coe woe aé [N(&,0,0) -| si veN(r')av], (134) 
O O 

Be = Pu f at [N(E,0,0) +f si v-N(r')av]. (135) 
O O 


Eqs. (134) and (135) give the general solution of the problem of amplitude and phase fluctu- 
ations for the case of large-scale fluctuations (ka >>1) under the additional condition that 
the distance is large compared to the scale of the inhomogeneities (L >>a). The form of the 
solution depends in an essential way on the size of the dimensionless parameter D = linjkee 
defined as the ratio of the size of the first Fresnel zone to the scale of the inhomogeneities. 
Following Gorelik Ba we Shall henceforth call this parameter the wave parameter. 

Eqs. (134) and (135) can be studied without specifying the form of the correlation coef- 
ficient N(r' ) in the following two cases which correspond to limiting values of the wave 


parameter: 1) D>>1, and 2) D<<1. This study reduces to an evaluation of the integral 


7) si v-N(r')av (136) 


in the two limiting cases mentioned. 


20. The Region of Large Values of the Wave Parameter (Fraunhofer Diffraction). The inte- 
grand in (136) is different from zero for values of r' which do not exceed the correlation 
distance a’ in order of magnitude. In this region of relevant values of r', the argument v 


of the integral sine does not exceed the value ia i.e. does not exceed the value 1/D 


x 
¢ 
oe 


Glnye 


Th 


In the case of a large wave parameter, v will be small. Therefore 
: 1 
Sl Vase-— , 
2 


Then Eq. (136) takes the form 


r=-% [ N(r')av. (138) 


The integrand N(r') has its maximum value N(E,0,0) for v = 0 (jn = ¢ = 0). From this maximum 


value N(r') falls off practically to zero in a distance vw 1/D, so that 
ee N(&,0.0) 
5 N(E,0.0). (139) 


The integral I can be neglected compared to N(£,0,0). Then Eqs. (134) and (135) can be writ- 


ten as one formula 


co 


Bo ge f N(E,0,0)dE (140) 


O 
or, if we go over to the dimensional variables L and &* = &/k 


oo 


Se eae feral f N(E*,0,0)de*. (140a) 


0 


Thus, in the case of a large wave parameter (D >>1) the mean square amplitude and phase 
fluctuations are the same and increase in proportion to the distance. For example, specifying 

ee 
e 


the correlation coefficient in the form N(é) = , we obtain 


uke. (141) 


(2 


Mintzer [21] , in studying pressure fluctuations of spherical wave sound pulses, used the 


method of small perturbations (in the usual form) and obtained the following formula: 


09 


Jax o 
—5 = ayok"L f N(E,0,0)aE. 
Ip | A 


The right hand side of this formula is twice as large as the right hand side of Eq. (140a). 
The reason for this difference is that the fluctuation of the field (pressure) is caused not 


only by amplitude fluctuations but also by phase fluctuations, i.e. 


2 
lAp | ne , =) 


=e 
Ip | 


In fact, the field can always be represented in the form 


_ _BtiStikx 


Expanding this expression in powers of the small quantities B and S and retaining only linear 


terms in the expansion, we find 


p =(1+B + is) en 
Averaging the field over an ensemble of media, we obtain 


ikx 
e 


for the mean field. Moreover, we have 


= : ikx 
PSPS PSR Ase ys 
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from which it follows that 


2 
| Ap | 3 7 
=B +5 
—je 
[>| 
ade Ne ; 
Since B =§ in the region of large values of the wave parameter, we have 
e 
| Ap | =) 
= 2B , 


1.e., the field fluctuations are twice as large as the amplitude (level) fluctuations. 
In a later paper Mintzer [22] found a criterion for the applicability of the first approxi - 
mation in the method of small perturbations, by requiring that the second approximation be 


small compared to the first approximation. This criterion is 
e 
uk aL <li, 


and agrees with the criterion indicated in Section 13 for the applicability of the scattering 
formula. This restriction does not apply to Eqs. (140) and (141), which were obtained by a 
different method. 

In studying the problem of the amplitude and phase fluctuations of a sound wave propagated 
in a turbulent atmosphere, Tatarski [32] used a scattering formula which is valid for large 
values of the wave parameter (Fraunhofer approximation). Therefore the formulas which he 
obtained for the mean square amplitude and phase fluctuations are suitable only in the case 


* 
of a large wave parameter . 


Ql. The Region of Small Values of the Wave Parameter. As explained in the preceding sec- 


tion, the region of relevant values of the variable v is of order 1/D. In the case of a small 
wave parameter (D << 1) the region of relevant values of the variable v extends to large values 


of v (v >>1). The problem reduces to evaluating the asymptotic value of the integral I in 


* In a later paper Tatarski [56] studied the same problem in a more general way. 


(7 


in the limiting case of large v. To do this we integrate (136) by parts 


ON(r') 


So dv (142) 


[6 @) 
(v si v + cos v) N(r') - f (v si v + cos v) 
O 


O 


= - N(E,0,0) - [ (v si v + cos v) ont") dv. 
O 


To estimate the asymptotic value of the expression v si v + cos v and certain other expres- 
sions, it is necessary to have an asymptotic expression for the integral sine. We can get 


such an expression from the asymptotic expression for the exponential integral Ei [39] 


k=n 
eee se". ChE! ae, (145) 
ea] x 
where 
n! d 
i 2 o 
IR | SS = & KS |x| er . d <1 4 
7 |x|" cos /2 
and from the formula 
Bi ei)" S161 eS) a eat ae, (143a) 


To evaluate x si x + cos x it is sufficient to take two terms of the series (143). Then we 


get 
xX ~ ae (144) 


1 ~ ———. . 
x si x + cos x = (145) 
The derivative ON(r')/dv is of order DN(€,0,0) and the integrand in (142) is of the order 
C 
D’N(E,0,0). Therefore, the integral in (142) can be neglected compared to N(E,0,0). Then 
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using (134) we obtain 


oe) 


=z 
S* = our" Fi N(E,0,O)dE, (146) 


O 
for the mean square phase fluctuation, or going over to the dimensional quantities L and &* 


co 


s-< aia ae f N(E*,0,0)ae*. (146a) 


O 


Eq. (146) differs from the analogous formula (140) for the case of a large wave parameter by 
the factor 2. The mean square phase fluctuation grows in proportion to the distance both in 
the case of a large wave parameter and in the case of a small wave parameter. 

Turning our attention now to the amplitude fluctuations, we note that in substituting 
(142) into (135) the leading terms N(§,0,0) and -N(—E,0,0) cancel each other. This compels us 


to consider terms of a higher order of smallness with respect to the parameter D. Then we 


obtain from (135) 


co co 


aoe ue f dé f (v si v + cos v) a dv. (147) 
O O 


Thus, in the case of a small wave parameter (D << 1) the amplitude fluctuations are an order 
of magnitude (with respect to the parameter D) smaller than the phase fluctuations. Since we 
wish to get a more suitable expression for the amplitude fluctuations, we integrate (147) by 
parts 


co 


ON(r' 
= [v(v si v + cos v) + sin v| ome") (148) 


fo si v + cos v) oN(z") dv = 
O 


O 


oo 
ent rt 
- f 3 DW stv + cos v) + sin SME) ay = 
Ov 
O 
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00 


ye 
( v [ov st v + 008 ») + ain] =F cos 2 


O 


+ 


r 5 
f € v[v(v si v + cos v) + sin v | - ; cos ) 3 dv 
O 


; -4 ae |, a jie v[v(v si v + cos v) + sin v| 


We can get an asymptotic expression for the curly brackets if we keep four terms in the 


expansion (143). Then we get 


de. we 1 6 
: sels ogi Ngee cans ee i dee ee ; a 
si v (3 =) cos v € <) Sin v, (149) 


2 v L(y si v + cos v) + sin v| - cos vw —— , (150) 


“N 


ee east is of order 


O/ 


The integral in (148) is of order D’N(E,0,0), and the term 


WIFE 


C 
DN(§,0,0). Keeping the leading term, we obtain 


(ee) 


; Pat 
f (v si v + cos v) ou(z’) dy = 5 oo 


0 


(15) 


2 
F : : C 2 ' 
Bearing in mind that v = ‘¢ + € )/4L", we can express | ae | in terms of deriva- 
3v (£,0,0) 


tives with respect to the variables n and €. In fact 
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Fy ek On Pe oN 
“EO Te 2.2 
where 
face, 
Tee ee 
om” oe 


Pen = (ES eS 


Therefore 


5 [ eae — 


Finally, Eq. (147) takes the following form: 


B= Lt i [PVN] eto dk. 
O 


In going over to dimensional quantities it preserves the same form, i.e. 
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(152) 


(155) 


(154) 


(155) 


(155a) 


Eqs. (146a) and (155a) coincide with Eqs. (69) and (81) of Part I, which were obtained 
using the ray representation. Therefore, the case of a small wave parameter corresponds to 
the geometrical approximation. Ellison [31] showed this in his study of the formula for the 
mean square intensity fluctuations, and Tatarski [33] also showed this in his study of the 
formula for the (unaveraged ) amplitude and phase fluctuations. In the geometrical approxima- 
tion, the mean square amplitude fluctuation grows in proportion to the cube of the distance. 

Tt is easy to give a qualitative explanation for the difference in the distance dependence 
of the mean square amplitude fluctuations in the regions of small and large values of the wave 
parameters. In the first case, the dimensions of the first Fresnel zone are small compared to 
the scale of the inhomogeneities. The deviations of the refractive index pw from its mean 
value within the zone all have the same sign. Therefore, all the waves scattered by the 
different elements of the first zone arrive at the observation point in phase, and the fluc- 
tuations increase rapidly with distance (like 1’). In the second case, the dimensions of the 
first Fresnel zone are large compared to the scale of the inhomogeneities. The deviations of 
the refractive index uw from its mean value have different signs at different points of the 
zone. Therefore, all of the elementary waves do not arrive in phase at the observation point. 


They partially interfere, and as a result the fluctuations grow more slowly with distance 


(like L). 


ee. The Region of Intermediate Values of the Wave Parameter. We can study the solutions 
(134) and (135) in the intermediate region by specifying the form of the correlation coeffi- 


cient. Let us set 


and introduce the variable v defined by (132). Then we obtain 


2 
rt 
N= oxp|- a = | : (158) 
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Substituting (158) in (134) and (135), we reduce the problem to evaluating the tabulated 


integrals [39] 


co co 


Cie 
| xG,0,0)a8 = { es /2 dé = — a‘, 
O O 
n 2 
; -iL' v/a! 
f Si ve v/a dv = peed ee arc tan Lt oa arc tan D 
bE" /at® gre > 
O 
where 
t 
a! ka 


is the wave parameter. Eqs. (134) and (135) become 


| 
ms 


unk aL (2 + = arc tan D), (159) 


. 
rls 


le eee - = arc tan D). (160) 


Obukhov [12] obtained these formulas in this form. In the limiting cases of large and small 
Wave parameters they go over into the corresponding formulas of Sections 21 and 22. 

All that we have said about fluctuations applies equally to electromagnetic waves. In 
scattering by large-scale inhomogeneities the important role is played by the region of space 
located in the immediate neighborhood of the straight line connecting the transmitter and the 
receiver. In this region the angle X between the direction of the electric vector and the 
scattering direction is close to a right angle and the factor sinx is close to unity, as a 
result of which the distinction between longitudinal and transverse waves disappears. The 
problem of fluctuations of the electromagnetic field caused by a turbulent atmosphere was 


investigated in the work of Rice [ 38] , Me gaw [18], Voge [20] and others. 
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Chapter VI 


CORRELATION OF FLUCTUATIONS 


The mean square amplitude and phase fluctuations still do not give a complete character- 
ization of the statistical properties of the wave field. The statistical properties of the 
fluctuations of the wave field can be characterized more completely by using correlation func- 
tions. However, the question of the correlation of the fluctuations of the basic character- 
istics of the wave field, of both theoretical and practical interest, has not received much 
attention. In this connection, three questions arise naturally: 

1. Is there correlation between the amplitude and phase fluctuations at the receiver? 

2. How does the autocorrelation between amplitude (or phase) fluctuations depend on the 
distance between receivers? 

3. What form has the time autocorrelation function for amplitude (or phase) fluctuations 
at the receiver? 

As far as we know, the first question has not been considered before; we answer it here. 
The second question has been examined by some investigators [6,14], but only in the ray 
approximation for the case where both receivers are situated in a plane perpendicular to the 
direction of wave propagation (transverse autocorrelation). Moreover, they have not dis- 
covered the simple law to the effect that the transverse autocorrelation between amplitude 
(or phase) fluctuations extends over approximately the same distance as the correlation 
between the random inhomogeneities of the medium itself, if the scale of the latter is large 
compared to the wavelength. Here we show that this simple law is valid both in the region of 
small values of the wave parameter (the ray approximation) and in the region of large ana 
intermediate values of the wave parameter. Moreover, we show that the longitudinal correla- 
tion extends over a distance which is considerably larger than the transverse correlation. In 
fact, if the distance between the receivers does not exceed the distance over which the ray 
approach is appropriate, then the amplitude (or phase) fluctuations are practically completely 
correlated. 


The third question concerning the time autocorrelation was investigated by Mintzer [23], 
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who assumed that the correlation coefficient of the refractive index fluctuations can be 
factored into two terms, one depending only on the coordinate differences and the other 


depending only on the time difference, i.e. that the correlation function has the form 
NCE, me So? Aa: a No? C. a Co) M(t, = ty)- 


This assumption is not justified in the case where the time change of the refractive index is 
caused by the motion of the inhomogeneities, which occurs in the presence of drift and con- 
vection. In this book we consider the question of how motion of the inhomogeneities affects 


the time correlation properties of the acoustic field at the receiving point. 


23. Correlation of the Amplitude and Phase Fluctuations at the Receiver. We begin by 
studying the cross correlation of the amplitude and phase fluctuations at the receiver. To do 
this we determine the form of the correlation function BS, using Eqs. (106) and (107). Multi- 


plying these equations and averaging them, we obtain 


je Ore + co 

Ra 2 ' t ’ ’ t ' ' ' t ! ' 

BS = u LLL $,(L'- 1,01) O(L'- 63,04) N(r')dgidésanidngdtiats. = (161) 
O O =O 


Considering only the case of a statistically isotropic medium, we shall suppose that the 
correlation coefficient N depends only on the magnitude r. Introducing relative and center of 


mass coordinates in (110) and (111), we obtain 


+ oo 


BB = LP Prey fv -e, Vaden? + Gea] x (262) 


e t t 
x 9, cs - §4) (A aa) 4 Gre |e" yas jas sanetava. 
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In this formula we can carry out the integration with respect to the variables y and z (see 


Appendix II): 


+ 
[| a,|u ee Ja + y)° + (2 + | x 
sar “tp Y@-v? +(3- fo 


-3 (0, a3 (65 + E4),0] + O(E1 = i420) , 


where p = Vine + (=, Then we obtain 


L'L'+ o 
Ra 12 t t t ! ! t 
BS = 5 Ul {Tf §,[2b" - (1 + &4),0] N(r')agjagsanae 
00-7 @ 
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Further simplification of Eq. (164) is possible in the case L > a. 


(163) 


(16+) 


Introducing the 


coordinates 2x = a + ES and &€ = = - ED » in this case we can integrate with respect to & 


between the limits - » and +o. Eq. (164) can be written in the form 


L! + © 

— 12 

BS = 5 [ ax inhi 9, [2L' - 2x, 0] N(r')dédnat 
O - © 


ED 


50 a ax ry $,(6,0) N(r')aganag. 
: a. 
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(165) 


Since N(r') is an even function of §, while $5850) is an odd function of §, the second inte- 
gral vanishes. In the first integral we can integrate with respect to x. Thus the problem 


reduced to calculating the integral 


i L' 


2 
i p 
rt oo = . 
f O,[2u" - 2x, e]ax f Bu(OL® = Ox) °° Bar — ox) 
O O 


: = fe) 
The substitution z = BOL? 2 Ox reduces it to the form 


(see the definition of the integral cosine [39]). Eq. (165) then takes the form 


BS = - = 2 fff ci (é) N(r')déandt - (166) 


Transforming to polar coordinates (p,@) in the (n,¢) plane, we obtain 


BS =- 


Mole 


Le oe) [e@) B 
uc f at [ ct GE N(r') pdp. (167) 
O Oo 


Finally, we introduce as before the variable v: 


O (168) 


Then 


B-- wu f dé f ci v N(r')dv. (169) 
O e) 
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Eq. (169) gives the general solution of the problem of the correlation of smplitude and phase 
fluctuations at the receiver. 


Suppose the correlation coefficient has the form 


2 2 
N(r') = exp[- pee Sy, (170) 


According to (158) the problem reduces to calculating the integral 


00 


f ci v exp “4 v | dve 
a! 


O 


We find in the tables [49] that 


CO 


2 2 
1 t hyt 
ff ctv em “Ss v J av = - By toe 1 + (5) ° 
a! a 


O 


Finally we obtain 


BS = aa uke” log (1 + D“). (171) 


In probability theory [44] the correlation coefficient is defined as the ratio of the 
correlation function of the fluctuations to their rem.s. values. Applied to our case this 


means that the correlation coefficient Rs of the amplitude and phase correlations has the 


form 


fee. He er (172) 


Using the formula (171) for the correlation function and formulas (159) and (160) for the mean 
Square values of the amplitude and phase fluctuations, we finally get the following expression 


for the correlation coefficient Ris? 


log (1 + D~) 


os (arc tan D) 


a 
mM 

ll 
Mlb 


° (173) 


At small distances (D << 1), when the ray approach is appropriate, Eq. (173) gives 


a 
te. 2 /2 oo O60 
At large distances (D >>1), Eq. (173) takes the form 


log D 


Bs Zi D 


i.e., the correlation coefficient falls off with distance and approaches zero. Thus, the 
correlation between amplitude and phase fluctuations which exists at small distances vanishes 


at large distances. 


24, Longitudinal Autocorrelation of the Amplitude (or Phase) Fluctuations. We now con- 
sider the question of the autocorrelation of the amplitude (or phase) fluctuations at differ- 
ent receiving points. Let us assume that the receivers are situated along the direction of 
propagation. Let the coordinates of the receivers be (L},0,0) and (L!,0,0). The amplitude 
and phase fluctuations are given by the formulas (106) and (107). In these formulas we have 
to put L'= L! for the first receiver and L'= L} for the second receiver. Therefore, the auto- 


ae 2 


correlation functions for the amplitude and phase fluctuations are given by the following 


formulas: 
LL, + © 
S_s 2 ' 1 t ! ! 
eee ITSLTS Oi (Ly bye Py) Obs - Eas pa) x (174) 
Oo O ~- © 


x N(x" dg jd sdnjangab sac. : 
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LIL' +0 


BE, CPOE Gyo] - HY, 0h) Bolus - bf, 0g) * oie 
OO = 


x N(r')déJdésdnidnsacsacs ; 


In these formulas the upper limit of the integration with respect to BD can also be set equal 


to Lj: In fact, because of the directional character of the scattering, the waves scattered 


by the layer bounded by the planes x = Ly and x = L, are incident on the second receiver but 


not on the first. Therefore, these waves can be neglected in the calculation of autocorrela- 


tion functions. 


setting the upper limit of the integration with respect to BD equal to Li, we obtain 


nieaae WP + © 


cle 
see PLL LL Bey - tp op ley - eg op » (276) 
oO O - © 


x 


N(x" Jag désanjansatiats , 


wee ve + © 


bagel 


bd 
bd 
II 


x 


N(x" ag désdnjangacrats : 


Using Eqs. (110) and (111) to transform coordinates and using Eqs. (114) ana (115) to carry 


out the integration, we find 
6s Ge (ie. or 
1°22" 2 i =r (176a) 


i 
BBL =54 (1, - 1), (177a) 
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where 


bi, 4.00 


al rey $3 - Ly 7: )s op] N ‘de dé sands, (178) 


O 0 = © 


KH 
lI 


LL, +00 


2 irre $, [ES + Ly (Es + BD ys 0 N(r’) dg 1g! sande. (179) 


O O - @© 


tH 
I 


These last formulas constitute a generalization of Eqs. (118) and (119) for the case Li} £1! 


The way in which they are manipulated further is analogous to the manipulation of Eqs. (118) 


and (119). Introducing the coordinates & = Es : Et and x = S(t + E1) and assuming L >> a, 


we obtain 


2 
t 1 > @) 1 
= Ly f dé 3 AL. & E eel ON a Mee Nir \pdp. (180) 
— @~ O 


Eq. (180) is analogous to (124). By introducing the new variable 
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Consequently, 
oO 
2 2 
I, = - f dé f si AEs ma - si it N(r' )pdo. (181) 
O O 


This formula is analogous to (125). 
Suppose that the correlation coefficient N(r') has the form (170).. Introducing the vari- 


able q = 0“ /2(ALt - E), we obtain 


1 = 0H 4 oan cs 
P wes ia aby oe |- EE oa oe 


3 
ce 
a 


a) » {are 


a! 


C 
2 e(AL' - €) exp |[- E /a'*] 
= exp|- >= sin q exp see x{dq = ————————~ , for AL' >€é. 
O 7 
For AL' < & we get the same result by setting q = ee ae E). By (180) we have 


*.° — expf- e°/at®] at 
T ian Tee f * (183) 


HE oO oC 
Sie. all Sai 2 eee) 


a! 


Calculation of an integral like qT, was encountered at the end of the preceding chapter (see 


p. 83). Repeating similar calculations, we obtain 


oe + L! 2 L! + t 
t = Ja at = ee 1 arc tan ares, 
2 2 OL? + L! D (18+) 
i. o a! 
gts 
1 CAL! 
- arc tan P 
2 Ly + Ly ate 
at= 
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receivers. Then 


tent 4 gy! 1 
a(Li + L3) ~ Li, 


and Eq. (184) takes the form 


7 1 i 1 AL! 
I, = Vn a'Li| = arc tan D - 5 are tan ( i > ; (135) 
where 
att 
Ses 
ie kaé 


In the region of large values of the wave parameter D the integral I, can be neglected com- 


2 
pared to the integral I: Then we obtain 
+ 2 
ere ee een i. 3 exp[- E /at*Jae a 
BB. = 6.5. == 1 Lh: ——_-———_—— ——_——_——_——_ . ab 
Le le 2 iL AL'<- OAL'e + Ee ( ) 
“op +4 Se Se 


a! 


If the distance AL' between the receivers is of the order of the correlation distance a', then 
the second term in the denominator can be neglected as compared with unity. In fact, the 
integrand differs from zero in a region of values of € of order a'. In this region of signi- 


ficant values of & the second term in the denominator will be of order 


a Z dle eo 1 
a! are (ieeye 


for ka >>1l. In this case (AL ~ a), we obtain 
LH 6k aL, 


ae ba aoe Vx 2.2 
2 


which coincides with the formula (141) for the mean square amplitude and phase fluctuation. 


Thus, at a distance AL of the order of the correlation distance a of the refractive index, 


2? 


both the amplitude and the phase fluctuations are practically completely correlated. 

We consider next, therefore, distances AL which are large compared to the correlation dis- 
tance a, i.e. AL >> at. In the region of significant values of € the terms AL'E and Bo have 
the orders of magnitude AL'a’ and Aca respectively. Therefore, we can neglect these terms 


as compared with the term (ant)*. Then Eq. (186) becomes 


“ te ea 
Bey = 684, = ere (187) 
tee. oe 
ka 
By (187) and (141) we obtain 
R= R, = ee (188) 
CAL . 
1 + 3 
ka 
for the correlation coefficients 
. ie alae 
iy = reer 9, ego ee 
BO g= 


As AL grows, the correlation coefficients fall off. Taking for a correlation distance the 
distance AL, | at which the correlation coefficient is reduced by a factor of two, we obtain the 


formula 


5 = (189) 


for the definition of the correlation distance. Comparing the left hand side of (189) with 
the expression for the wave parameter D = ieee we arrive at the following conclusion: 


The longitudinal autocorrelation of amplitude (or phase) extends over the distance within 


—t ee ee OO 
anwRPNl eeu —_SEeee ee eee eee eee 
=a EER Cee ee eee —ee 0 Ee ee eee eee 


—eeieD eee ee 


Because of what we have just shown, it is no longer necessary to consider the other limit- 
ing case, corresponding to small values of the wave parameter: in this region there is com- 
plete longitudinal autocorrelation of the amplitude and phase fluctuations. Moreover, it fol- 
lows from the preceding section that in this region the amplitude and phase fluctuations are 


orn 


completely cross-related. 


25. Transverse Autocorrelation of the Amplitude and Phase Fluctuations. Suppose that 


both receivers lie in the plane x = L at a distance 2 from one another. Let the coordinates 
of the receivers be (L,0,0) and (L,0,2), respectively. The amplitude and phase fluctuations 
are given by Eqs. (106) and (107) both for the first receiver and for the second receiver. 


The difference is only that we must take P} = Jn + a instead of o' for the first receiver 


and PS = V15 + (vt! - Co) for the second receiver. Therefore, the autocorrelation functions 


for the amplitude and phase fluctuations are given by the following formulas: 


Dee + © 


meee [PLP ar - a op Bar - ty, op) x (190) 
00 -o 


N(r') alae sdnjdnjats ass, 
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hae + 0 


BB, ou FLL GY Q(L' - 84, ef) O(L" - &5, 04) x (191) 
Oo O - © 
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N(r' )déjagtdnjdnjacidcs. 


Using Eqs. (110) and (aaa) to go over to relative and center of mass coordinates, we obtain 


be + co 


mB [fee Mev Et], om 
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N(r' )agjaésdndgdyaz, 
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Tea + 0 


a Ffiiiab. reo 


<5» ae (A ~ y)© + (4! +2- | x 


x N(x! )désd&Sdndg dydz. 


The integration with respect to the variables y and z can be carried out in these formulas 


(see Appendix II): 
+ 7 
[f3fu wy Veo eee P| aoe 
70 oo a Vi (2 = y)* + (2! +2 = oP av 


= (5, (83 = ae p) + , [2u' 5 (E+ + Es), 3) J 


+ © 
[| g,| oe (A Es y)* + (2 y | x (195) 


t t e t 
<3, - £4, (A -y) + (2' + 2 7 oP at 
L 
- 5 (4 (64 ae oe ee ea é}) , 
where i = a + (C+ Bey. Then we obtain 


I, + Ly), (196) 
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—— 1 2 
fia 2 © ea aa) (197) 
where 
L'L! +0 
I, = LTT oe) = 5» P) N(r' )dandtagsdes, (198) 
O O = 
Li! + oo 


rH 
II 


2 LTT 9, [eu' - ae) N(x! )dnagag tags. (199) 
OO -—o 


Eqs- (196) - (199) give in general form the solution to the problem of the amplitude (and 
phase) autocorrelation at different points of the plane x = L. They are identical in their 
structure to Eqs. (116) - (119) of the preceding chapter, which give the mean square amplitude 
and phase fluctuations at the receiver. The only difference is in the definition of pO. Pre- 
viously, p was defined by the formula Fe = 1+ ae but now po is defined by the formula 

ae = noe (C + guy: This seemingly small difference leads to a destruction of axial symmetry, 
which makes it inappropriate to introduce polar coordinates in the (n,¢) plane and makes it 


impossible to make a general study of the integrals (198) and (199) in the way that we 


studied the integrals (118) and (119) for the limiting cases of large and small values of the 


wave parameter. Instead we have to start with some explicit correlation coefficient N(r') and 
carry out all the calculations in a somewhat different order than before. It can be assumed 
that the qualitative conclusions which are obtained for a special form of the correlation 
coefficient remain valid for another choice of the correlation coefficient. 

If the distance is large compared to the correlation distance (L >> a), then introducing 
the relative coordinate —& = EF - ES and the center of mass coordinate x = =(e! + Es), we can 


integrate with respect to & between the limits - ~ and +o. Then we obtain 


00 


Ae + 00 


I, = [ dx f de fi $, (£50) N(r')dnat, (200) 


_ = © 
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Oe f dx f dé {| g, (aL! - 2x, po) N(r')dnde. (201) 


Suppose the correlation coefficient has the form 


2 2 12 
N(r') = oc] St), (202) 


a 


By (101) the integrals I, and I, can be written in the following way: 


7 r yp al 1+ (C + ut)e Em 4 n+ ge 
I, = f ax f dé {| nk sin BE exp |- 2 dndc, (203) 

O - © — 0 a 

Ls + 00 + 0 

O o 
7 ale : + (C + Lt) ‘ 

I, = [ dx [ dé {| En(L'- x) sin ee x (204) 

O - © - © 

at ey r= 
x exp 2 dnag. 
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Consider now the calculation of the integral I,- First of all we carry out the integration 


with respect to the variables nf and ¢: 


+ © 
. o o O o 
gt 
ff sin THe tO) exe At dna (205) 
a a! 
a 2 12 2 3 
nial, ne + (6 + t') 
bes a. io 2 2 
aS eG 2 ke eG a) 
A {I exp ae a BE danas 


98 


+ co 


+ 00 
e O Cc 
-i i aN one C (C+ tt 


- © 


+ the complex conjugate integral. 


Moreover, 


In order to integrate with respect to €, we introduce a new variable €' by the substitution 


alle 


S Sh) oe or (C' =€ +c - id, 
=~ -i 
se 
a 
where 
OLie 
Lt 7c 
Cae 4 Ate 5 (206) 
1 + us 1 + oA: 
we a 
a a 


The quantities c and d are finite for any value of &. The new variable €' varies from 


-a.- id to +» - id. The integral with respect to € can be transformed into the form 


+o-id 
12 fe ° 
exp , tle [ exp ‘(3 - 25) ptr d¢'. (207) 
(2¢/a'”)-i ee a" 


The integration along the line parallel to the imaginary axis (a distance -id away from it) 


can be replaced by an integration along the real axis. Indeed, choose the contour of integra- 


tion as shown in Fig. 4. Then by Cauchy's theorem 


fefefefee ne 


I II Iii IV 
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Fig. 4 The contour of integration for Eq. (208). 


On the segments II and IV, where (¢' Ste yy We have 


ha 4 
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xy 2 ate - L x. 
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For sufficiently large XO the term x y/6 can be neglected as compared with the term xo fa’ 
O 


i.e. 
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Therefore > 0 as xX, @, and from (208) we obtain in the limit 
rr, ty 


I LIL 


Integrating along the real axis in (207) we obtain! 


are/ ar 
(2 /a'®) - 4 


Ae 


exp 


ste 
ee ee 


By using (203) the integral I, can be written in the following way 


1 
L' + co y yx 
= ee a e I s J. 
ae ne oe es (203a) 
A vee (2e/a'“)-i (2#/at)4i 


where 


" a 
ADRS?) ot 


O24 is the complex conjugate of 7): In just the same way we carry out the integration with 


respect to n and € in (204), obtaining 


L! +00 y y% 
Eat = f ax f ae of /at ee es 5 4 , (209) 
: O ae [4(L'-x) /a' ‘|- 1 [u(Lt-x)/at“]+ i: 


T The result can be obtained more quickly if in doing the integration in (205) we use the 
formula 


+o 
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where 


ino ne 
ivwtSa (210) 


"2° fh (Lt-x)/a'<]-i 


x 
(%, is the complex conjugate of Yo) 


We continue the calculation of the integral I_,. Setting a A + iy, » we find 


1 


Gt, ic Sar vit = (2! /a"©) (28 /at®) 


y= 
1 2 1 z 
1+, 1+, 

a a 


+o cos yi a sin "4 
E a! 
= t ee t 
I, =1L f exp [. - al exp [7+] a dé. (211) 
= 1+—7 
co ne 


The first factor in the integrand is appreciable only for values of & which do not exceed a! 
in order of magnitude. In this region the term at /at® can be neglected as compared with unity, 


since a' >>1. Thus, assuming that a' >>1, the integral (211) becomes 


ce ie 2 
L 2é E 
= 1 eae cence = 
I, =. exp [- = f so ( 3 3] exw 3 dé , (212) 
a oe a a a 
and finally 
2 2 2 
2 L L 1 
I, = x Lak” exp[- al exp (2 7 ‘ (213) 
a a 
As a result of the change of variable y = k(L'-x)/at®, the integral I, takes the form 
D 2; 2 ge ge : 
: Va 3 / exp[(it"/a")/(y-i)]  exp[-(i2°/0°)/(y+t J] 
yg ee a se ti 214 
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; ? . 92/2 : 
If we introduce the new variable z = (it"/a~)/{y-i), then the first integral in (214) becomes 


eye 
P exp [(i0?/e”) /(y-2)] “oe ok 
f you dy = == GzZ% 
0 (10° /a°J/(D-1) 


; oie 
After the change of variable z = (-it Jay / GRD the second integral in (214) becomes 


Wenig rag ace 
» exp(-(it*/e")/(y+i }] 2 
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° 17 /a® 


Therefore the integral I, can be represented in the following form: 


2 
Vit , 
I, = = a'?[Bi(e,) ‘ Bi(e,)], (215) 
where 
O42 2; 2 
—  aith/e _ ita 
2 Spe eS ST “pe (216) 


and Ei denotes the integral exponential function. 
First of all we shall show that our formulas give the familiar expressions for the mean 
Square amplitude and phase fluctuations in the special case 4 = 0. In this case, we find from 


(213) that the integral I, is 


ine Sak ay. (217) 


Turning to the integral I,; we use the following series representation of the integral expo- 


nential function Ei [39]: 


co 


k= . 
Bi(x) = C + log x + ; 
( ) g k- kl 
=) 


k 
Then we obtain 
j co K co eK 
Ei(e,) - Ei(e,) = log —_ an : 
o 1 aT k-k] kk] 
k= k=1 


For £ = 0 all terms of the series vanish and we obtain 


- a. 
heen [Ei(e,) ~ Bi(e,)] = 1o8|- 7 ql = 2i are tan D, 


and 


IT, = va Ka? arc tan D. (218) 


Finally, using (217) and (218), we obtain the familiar formulas 
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UL %“al.(1 += are tan D), (219) 
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Returning to the question of the transverse autocorrelation of the amplitude (and phase) 
fluctuations, we investigate the corresponding amplitude autocorrelation coefficient R. and 


phase autocorrelation coefficient Ro: beginning with their definitions: 


a= 5 ; (221) 
=p 
ms 

ae ets (222) 
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Using Eqs. (196), (197), (213) and (215), we obtain 


2 2 2 
exp[- 4] exp[- € =| ] - spFile,) - Bi(e,)] 


Re, (223) 


1 - = are tan D 


o 2 2 
exp |[- =] exp |- (= = 7 + ss5[Bi(e,) ~ Ei(e, )] 
a a 
R, (£0 —QvQV$V\Q"®\$") ————— (224) 


1 += are ten D 
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As is apparent from these formulas, the autocorrelation coefficients depend only on the three 
dimensionless parameters t/a, D, ka,where the dependence on the parameter ka can be neglected 


(in the case being considered ka >> 1) by setting 


c ka os a 
a 


Eqs- (223) and (224) then become 


R, = 7 ; (225) 


Ll - D arc tan D 


eo [=5] + arplEt(e,) - Ei(e,)} 


L +e are tan D 


(226) 


In the case of small values of the wave parameter (D << 1), these formulas can be simpli- 


fied. In fact, from (216) we get the equality 


2 
s ot 
-2i 2 D 
€ ae ; 
= F 1 + D° 
from which we see that the difference €, - €, is small provided that ai ey << 1. The last 


2 1 
inequality will be met, for example, under the condition that D << land twa. If it is met, 
then we can confine ourselves to the leading terms in the expansion of Bi(e,) - Ei(e,) in 
powers of small increments of the argument E., - €,- The expansion procedure is simplified if 


we introduce a third point aE = ia oe lying halfway between the points ey and En and if we 


represent the difference of the integral exponential functions in the form 
[Ei(e,) : Ei(e )] - [Ei(e,) - Ei(e )]. 


The quantity in the second set of square brackets is the complex conjugate of the quantity in 


the first. Therefore, if we expand one of these quantities in series, we automatically 
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obtain the expansion of the other at the same time. We expand the first bracketed term in 


series about the point Ey nee = 


ade east 2 Te Gopaitit > 
I = l = oe oe 
Ei(e,) Ei(e.) Bi Ce.) Ae + 35 Bi (e0) i Z Bi (e.) fem + 


Here 


fo) € 
; € ott HE ug O 
Ei'(e |) =——, EBi'(e_) = (2 haat 


€ 
ne 7 1 fa 2 O 
Ei (e.) = (2. ror 2.) e-, 
O O 


Oo € E 
2.2, 
aC ° ae 2° 2 iS 
See ee. Sp ae e (iD - D- iD’). 
? 1+ D a 


The expression for Ace is given up to terms of the third power in D. With the same degree of 


accuracy, we obtain 


€ 
: O|. aa a. pees 
Ei(e,) Ei(e)) e [ip - me + é,)D - Ae a ee cv), 


whence 


<6 D a Se 
Ei(e,) - Ei(e,) = 2iDe F - z(1 thE it ar 


In calculating the amplitude correlation coefficient R, using iq. (225), it is necessary to 
consider the quadratic term in D, since the leading terms cancel each other. In calculating 
the phase correlation coefficient R_ using Eq. (226), the quadratic term in D can be neglected. 
Correspondingly, writing 


2 


1 
l- gare tenD~> and 1 += arc tan DA2, 


we obtain 
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(227) 
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R, = exp(—>) - (227a) 
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The phase autocorrelation coefficient R. has a Gaussian shape and coincides with the 
correlation coefficient for the refractive index fluctuations. The amplitude autocorrelation 
coefficient R. is not Gaussian. However, in this case also, as can be seen from Fig. 5, the 
correlation between the amplitude fluctuations extends over a distance of the same order as 


the correlation between the fluctuations of the refractive index. 


Fig. 5 The transverse autocorrelation for the case of small D Go canis 
1 - the autocorrelation coefficient of the amplitude fluctuations; 
2 - the autocorrelation coefficient of the phase fluctuations, 
identical with the autocorrelation coefficient of the refractive 
index fluctuations. 
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In the region of large values of the wave parameter (D >> 1) 
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Using the formula 
Ei(-ix) - Ei(ix) = 21 5 - 2iSix, 


we transform the autocorrelation coefficients (225) and (226) into the form 
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The dependence of the autocorrelation coefficients on the distance 4 between the receivers is 
shown graphically in Fig. 6 for the case D = 10. The middle curve shows the correlation 
coefficient exp(- t/a) of the refractive index. From the curves it is clear that the trans- 
verse autocorrelation of amplitude and phase fluctuations extends to a distance of the order 
of the correlation distance of the inhomogeneities in the medium. We came to the same con- 
clusion in considering the case of small values of the wave parameter. Therefore, we can 
assume that this result also remains true in the region of intermediate values of the wave 
parameter (D~ 1). 

Of course, it is more difficult to investigate the dependence of the autocorrelation 
coefficients Ry and Ro on the parameter t/a in the intermediate region. The difficulty is 
connected with the fact that the integral exponential function is not tabulated for complex 
values of the argument. However, we can obtain the asymptotic values of the autocorrelation 
coefficients for sufficiently large values of the ratio t/a if we use the asymptotic expres- 


sion of the function Ei encountered in the preceding chapter. In fact, in this case 
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in absolute value. 


Sioa a D D aD D 2. COS a? {©XP 
a LD a a 


1 
a 
2 
a 1 v2 
2 


After some elementary manipulation, we obtain 


14+D L a 1+D 


ng -= arc tan D 
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exp(- a ee sin(—> 5) ere cos( z 5) exp/( 5 
1+D a 1+D 1+D 
Rn ee (251) 


1 + = arc tan D 


for the autocorrelation coefficients. Eqs. (230) and (231) contain only squares of the ratio 
t/a, which makes them applicable for t/a as small as 4. For such a value of the ratio t/a 
and for D~ 1, the autocorrelation coefficients are small as compared with unity. Consequent- 
ly, in the intermediate range, the autocorrelation of amplitude (and phase) fluctuations also 
extends over a distance of the order of the correlation distance of the inhomogeneities in 


the medium. 


Fig. 6 The transverse autocorrelation for the case of large D (D = 10): 
1 - the autocorrelation coefficient of the amplitude fluctuations; 
< - the autocorrelation coefficient of the refractive index fluc- 
tuations; 43 - the autocorrelation coefficient of the phase fluc- 
tuations. 
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26. The Quasi-Static Condition. Until now we have regarded the distribution of inhomo- 


geneities as static, neglecting their change as a result of heat conduction, diffusion, con- 


vection and drift. We can neglect this change if the propagation time t = L/e is small 


compared to the characteristic time scale of change in the inhomogeneities. However, if this 


condition is not met, then in calculating the amplitude (or phase) fluctuations at the time t, 
we have to take into account the refractive index fluctuations at the time t'=t - r/e, where 
ris the distance from the scattering element to the observation point. 


In this case, the 


basic formulas (106) and (107) are written as follows: 


L' + w 


[ff &e- E'y,o') wCE', nt, 6', "dé tant att, (232) 


Oon- @ 


5 L'3030,,5) 


alee + co 


[| $,(u'- ep) Ce ey ye Tye de dy ae! (233) 


QO - a 


B(L',0,0,t) 


where t' = t - r/c. Squaring and averaging (232) and (233), we obtain 
Leo + co 


aw PPL LL tar ap op Rar ey op) x (234) 
ome) - © 


x M(EL - Ebon} - nfoSf = Gbeth - eas jasg,anjangatiass 


—_ eee L'L' + oo 
cae PLL LL G,(L' - &], ef) O(L'- 63, p35) x (235) 
OO - © 
x (EI - Bb nt - aBsbr - b4,t' - t")dbsabsanjanzatidss,, 
where 


tt 
N(E! - §3, nL ~ ne $4 - $5, t' - t") = 


NGS ) : 


iil 


Here we assume that the correlation coefficient N depends only on the coordinate and time 
differences, i.e., we assume that the process is stationary in time and homogeneous in space. 


The time difference is 
tt at eS (es r,/c) -(t - eal) = (x, - r,)/¢- (236) 


The correlation coefficient N is different from zero if ry - Py does not exceed the 
correlation distance a in order of magnitude. This means that the time difference t' - t" 


does not exceed the quantity a/c in order of magnitude, i.e. 
et see fo ee a/c ° (237) 


The quantity a/c defines the time it takes the wave to go a distance equal to the correlation 
distance (the scale of the inhomogeneities). If this time is small compared to the correla- 
tion time T of the refractive index, then the time difference t' ~- t" can be set equal to zero, 
and any explicit time dependence in our formulas vanishes. Consequently, the quasi-static 


condition takes the form of the inequality 
a/e << T. (238) 


The quasi-static condition (238) for averaged quantities is very much less stringent than the 


quasi-static condition for unaveraged quantities, which clearly can be written as 
L/e << T. 


In all actual cases, the quasi-static condition (238) for averaged quantities is evidently met 
with a big margin. This justifies the quasi-static assumption which we have made from the 


very beginning. 


ef. The Time Autocorrelation of the Amplitude Fluctuations. The change of the inhomo- 


geneities in time produces a change in the frequency of the scattered waves and broadens the 
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frequency bandwidth of the incident radiation. The nature of this broadening can be inferred 
from the form of the time autocorrelation functions of the amplitude and phase fluctuations. 
Thus, for example, the frequency bandwidth is determined by the correlation time of the radia- 
tion. However, we shall not be concermed with the frequency structure of the fluctuations, 
and only the autocorrelation functions are examined below. 

An appreciable weakening of the autocorrelation of the amplitude fluctuations can be 
expected after a time interval Tt commensurate with the correlation time T of the refractive 
index. If the quasi-static condition (238) is met, then the time interval t is also large 


compared to the time a/c, Lee. 


t >> a/e. (239) 


We write Eq. (233) for the times t, and t,, assuming that they are separated by the interval T: 


i 2” 
L' +0 
= to gt ' t t t t t t t Oho 
B(Lt,0,0,t,) = f ff @(L' - €', of) w(et, nt, St, tt)agtantas’, (240) 
oO - @ 
L' +o 
t = = t t t : t " t tq er ayia 
B(i',0,0,t,) = fff BG" - £4, ot) wet, nt, St, tag tantat (2h) 
Oo - © 
where 
thetp-r/e, tt ety-2/c. (22) 


Denoting the time autocorrelation function of the amplitude by F(t), by definition we have 
F(t) = B(L',0,0,t,) B(L' ,0,0,t,,) : 


MeLtiplying (240) and (241) and averaging, we obtain 


1;'- 1;! + 00 
2 t t t oh 
a =e ff PP ff Bar - ep ep Be - ey 05) * (25) 
O O - © 
x N(ES ~ D? 7h = No? Ci ~ >? GS t")dg idk Sans dngde sas. : 
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By (242) we have 


ee ee gee (ce r,/¢) - (t, - r/c) (244 ) 


since (r, - r,)/e~ a/c. Therefore Eq. (243) becomes 


R(t) = uP Prrrry G(L' - et, pt) G(L" - eb, ph) x (25) 
Oo O ~ © 


t t t t t t t 
x N(r ,T)GE Fae Sanjantacsacs : 


Mintzer [23] started with the assumption that the correlation coefficient N(r',t) sepa- 
rates into two factors, one of which depends only on the coordinates and the other only on the 


time, i.e., that the correlation coefficient can be represented in the form 
N(r') M(t). (246) 


Then (245) becomes 


LL! + 0 


exer) Of PPE oat - ey, 0) Bye - eg 0g) x (247) 
oO O ~ 


x N(r' ag jdéSanjangacsags i 


Denoting the amplitude autocorrelation coefficient by R(t), we have 


R(t) = oo (248) 


B 


by definition. Then by (247) and (109) we get the very simple result 


R(t) = M(t) . (249) 
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Thus R(t), the time autocorrelation coefficient of the amplitude, coincides with M(t), the 
time autocorrelation coefficient of the refractive index. It is possible that Mintzer's 
assumption that the time and space coordinates can be separated is justified in the absence 

of any motion of the inhomogeneities which are produced by drift and convection. In this case, 
change of the inhomogeneities in time might be produced, for example, by turbulence, heat con- 
duction and diffusion. However, such a separation cannot be valid when drift or convection is 
present. In this regard, it is of interest to examine the case where the change in the 


inhomogeneities in time is caused by their motion. 


We assume that all the inhomogeneities move with the same velocity v, as a result of 
drift, say, and we assume that the change in the inhomogeneities results exclusively from the 
drift, while other factors (turbulence, heat conduction, diffusion) play no important role, 
i1.e., changes produced by these factors proceed much more slowly. Then in the coordinate 
system moving with the flow, the correlation coefficient depends only on the coordinates, i.e., 


has the form 


N(x, A Rp a AAP ey. = Zr) « 


Using the Galilean transformation formulas 


: . : 250 
Y 7 %¥o = 1 7 Nee (250) 


we obtain 


~£,- - -t,- 251) 
N(E, bo Vege 3 unl No? C1 Co vt) ( > 
in the coordinate system &,n,¢ fixed with respect to the receiver. (The coordinate systems 
are oriented in such a way that the flow velocity v lies in the planes xoz and tof.) It can 
be seen from (251) that in the case under consideration there is no separation of the space 


and time variables. 
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It is easy to see that in the case of a homogeneous flow the problem of time correlation 
reduces to the problem of space correlation which has already been solved. To see this it 


suffices to use the principle of relativity to go over to the coordinate system x,y,zZ moving 


with the flow. The receiver moves with the velocity -V with respect to this system. We assume 
that at time t the receiver is located at the point A and that at time t + T it is located at 

the point E, where the receiver displacement AE = ¢ satisfies the condition = VT. We denote 
the amplitude fluctuation at the point A at the time t by B(A,t) and the amplitude fluctuation 
at the point E at the time t + tT by B(E,t + T). Since the distribution of inhomogeneities can 
be considered static in the coordinate system moving with the flow, inasmuch as we are dealing 


with time intervals of order tT, the amplitude at any point does not change with time, i.e. 


B(E,t +7) = BCE, t). 
Therefore 


BUA, BCE, b "2 ). = BOA ty BCE yt): 2 (252) 


Here the left hand side represents the time correlation function F(t) in the coordinate system 
fixed with respect to the observer, while the right hand side represents the space correlation 


function P(t). Thus we have 
F(t) = Fj (2). (253) 


Since the longitudinal displacement in a time tT of order a/v cannot produce an appreciable 
decorrelation, in Eq. (253) it is enough to consider only the transverse displacement %’ =vit. 
Z, Z, 


Then we get 
(a) eRe): (254) 


We note that the correlation functions of the phase fluctuations also Satisfy the relation 
(25 yi. Tt only remains to use the corresponding results of Section 26, where the transverse 


displacement was denoted by £. Thus, for example, Eqs. (225) and (226) give 
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2 
ae exp |- 3 rd - 5 Bile.) - Ei(e,)] 
R, = (255) 


] - - arc tan D 


ve 2 i: 
exp|- 2 Tt] + Spee) : Ei(e, )] 
R(t) = aa bias (256) 
1 += arc tan D 
D 
where 2 
al “2 r 
= * 
Ee, =- -———, Ey Bye. s 
S fae 1 9. (257) 


In the Fraunhofer zone (D >> 1), they take the particularly simple form 


2 
R(t) = R(t) ~ exp oe “|. (258) 
Defining 
te = ) (259) 
we get 
e 
R(t) = R(t) = exp(- za (260) 


It is clear from (259) that the correlation time T depends on the transverse component of the 
flow velocity or (what amounts to the same thing) on the transverse component of the receiver's 


velocity. 


28. Comparison with Experiment. This section is devoted to discussing some acoustical 
and optical phenomena from the standpoint of the fluctuation theory presented in the two pre- 
ceding chapters and then comparing the theory with experiment. The first experimental studies 
of amplitude and phase fluctuations when sound is propagated in the atmosphere are due to 


Krasilnikov [52,53]. Later, using an improved method, Krasilnikov and Ivanov-Shyts [7,8] 
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also made experimental studies of amplitude and phase fluctuations when sound waves are pro- 
pagated in the atmosphere, and they compared the experimental results with the theoretical 
results obtained by Krasilnikov [6] using the ray approach. They discovered that the theore- 
tical result concerning the phase fluctuations was in satisfactory agreement with experiment, 
while the theoretical formula for the amplitude fluctuations was not supported by the experi- 
mental data. Fig. { shows a logarithmic plot of the dependence of the mean square phase 
fluctuations (over various time intervals) on the distance between the transmitter and 


receiver. 


~ Uv 
1.2 cS 
& 
o" s 
4.0 
0.9 39 a5 67 m 


Fig. ( Dependence of ra on L for the frequency v = 3 Keps 
(logarithmic scales): A,At = 0.04 sec; 0,At = 0.08 
sec; O ,At = 0.2 sec. 

(After Krasilnikov and Ivanov-Shyts) 


The line V2 is drawn through the experimental point corresponding to the distance L = 22 nm, 
: ; wae, We 1/2 

It is clear from the figure that the theoretical law, according to which VS ~L » agrees 

satisfactorily with the experimental data. The dependence of the mean square amplitude fluc- 

tuation V Bo on the distance L is shown in Fig. 8. 


The location of the experimental points with respect to the curves pi/2 5 13/2 gives an 


, . 2 
idea of the behavior of B™ as the distance changes. As can be seen from the figure, the 


amplitude fluctuations, like the phase fluctuations, grow with distance in proportion to pi/2 
whereas the ray theory leads to the stronger dependence Je ual of the amplitude fluctu- 
ations on distance. Trying to explain the discrepancy between theory and experiment, Krastl- 
nikov [7] (and others ) proposed that the ray approach is not suitable under the conditions of 


the experiment. Apparently, the cause for the discrepancy was correctly ascribed. In Lact. 
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} 


the condition that the wavelength be small compared to the scale of the inhomogeneities is 
not met: The inner dimension of the turbulence in the layer of the atmosphere near the earth 
is of the order of 1 cm, which is much smaller than the wavelength of the wave used in the 
experiment, namely 6.6 cm (5000 Keps). On the other hand, the diffraction theory gives the 


dependence of the amplitude fluctuations on distance which is experimentally observed, i.e. 
Veen me 
vp? 


O7 “ 
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0.5 ® 
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Fig. 8 Dependence of V Bo on L: 0, measurements of Aug. 1, 

frequency v = 5 Keps; A, July 31, v = 5 Keps; 

+, Aug. 51, v = 3 Keps; e, Aug. 9, v = 3 Keps; 

O , Aug. 10, v = 43 Keps. (After Krasilnikov and 

Ivanov-Shyts ) 

Sheehy [15] studied the fluctuations of signals propagated in deep water. In his experi- 

ments a series of sound pulses was sent at each fixed distance and the size of the relative 
fluctuations of the pressure amplitude was determined and expressed as a percentage of the 


mean amplitude (see Fig. 9). The experimental line (1 in Fig. 9) shows that the size of the 


fluctuations grows as the square root of the distance, although the scatter of the experimen- 


y 


tal points is large. Using Liebermann's data (ue =5x 10 
2; 2 
r /a 


, a = 60 cm) and a correlation 


coefficient of the form e- 
V= V eh! ? ‘e kal. 


between the relative pressure fluctuation V and the distance L (line 4) and also the depen- 


2 


» Mintzer [21] represented on the same graph the dependence 


dence V~ 13/ which follows from the ray theory (line 2). As can be seen from the graph, 


the latter dependence is in strong contradiction with experiment. This is explained by the 


119 


fact that Sheehy's data pertain to the Fraunhofer zone. In fact, for distances in the range 


50 m to 2500 m for a frequency v = 24 Keps and for a scale of inhomogeneities a = 60 cm, the 
wave parameter varies from 6 to 300, i.e., is much larger than 1. The results of the diffrac- 


tion theory (line 3) are in good agreement with the experimental data. We note that between 


; 2 : 
the dashed boundary lines, containing 90 percent of the data, the quantity ua varies from 


- - -7 
1.5 X 10 6 to 1.5 X 10 io (The experimental line 1 corresponds to the value 5 X 10 .) 


0 90 100 250 500 1000 2900 L,m 


Fig. 9 Dependence of the size of the relative fluctuation 
of the pressure amplitude of the direct signal on 
distance. (After Sheehy and Mintzer) 


In optics, because of the smallness of the wavelength, the wave parameter D = utes can 
stay small compared to unity even over large distances; as a result the region of applica- 
bility of the ray theory is greatly extended. Convincing support for this can be found in the 
phenomenon of the twinkling of stars. Butler [16] studied experimentally the brightness 
fluctuations of stars as a function of the zenith distance. The results of his observations 
are shown by the points in Fig. 10. The path of the ray in the atmosphere is L =h sec Zs 


where z is the zenith distance and h is the height of the layer containing the inhomogeneities. 
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Therefore in the graph the relative brightness fluctuation is given as a function of sec Ze 


Fig. 10 Dependence of the relative brightness fluctuation of 
stars on the zenith angle. (After Butler and Megaw) 


Megaw [17] compared these data with theory. As can be seen from the graph, the dependence 

(AL/1)* ~ (sec 2)°!* nolds for zenith distances which do not exceed 60°. At greater zenith 
distances the brightness variation obeys the law V (AL/1)* ~(sec 21/2, The graph is inter- 
esting in that it encompasses both laws at once and shows quite clearly the transition from 
one law to the other. Van Isacher [19] draws similar conclusions from Butler's observations. 

Experimental material on the correlation of fluctuations is extraordinarily skimpy, and 
that which is available requires correct interpretation. Let us examine some experimental 
hydroacoustic data pertaining to time correlation of amplitude fluctuations. Such data 
(obtained by Sheehy) is presented in Mintzer's paper [23]. Analyzing this data, Mintzer came 
to the conclusion that the time correlation coefficient does not depend on the distance. On 
the other hand, experiments carried out at the same distance but at different times give 
widely differing curves for the correlation coefficient. Since the experiments were carried 
out on moving ships, it is natural to assume that the difference in the curves was produced 


1e1 


by the difference in the ships' motion. Starting with this assumption, we can interpret the 
results of the experiments, using the theory developed in the preceding section. 

The experiment consisted of 50 series of pulses. Each series consisted of 50 pulses at a 
frequency of 24 Keps. In 16 series the length of the pulses was 10 msec, and in the remaining 
series it was 430 msec. The transmitter and receiver were at a depth of 15 m. Since the 
pulses are separated by one second intervals, the values of the correlation coefficient were 
found for T = 1,2,3,... sec. Fig. 11 shows the time 769) of the first zero of the correlation 
function as a function of the distance between the transmitter and the receiver. 
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Fig. 11 Time 760) of the first zero of the correlation function as a function 
of the distance. The circles correspond to pulses of length 10 msec 
points correspond to pulses of length 40 msec. 
(After Sheehy and Mintzer) 
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The points and the circles refer to the pulses of lengths 350 msec and 10 msec respectively. 

It appears that there is no particular dependence of 760) on distance. On the other hand, the 
variability of 760) for different series at the same distance from transmitter to receiver is 
striking. Since there is no dependence on distance, the correlation function can be found by 
averaging over all 50 series. The values of the correlation coefficient obtained as a result 


of this averaging are denoted in Fig. 12 by circles. 
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Fig. 12 Dependence of the average correlation coefficient 
on time. The experimental points were obtained by 
averaging over 50 series of sound pulses. 

(After Sheehy and Mintzer) 


The variability of the correlation coefficient can be explained by its strong dependence 
on the velocity of motion of the receiver, which is apparent from Eq. (258). The time depen- 
dence of the correlation coefficient shown in Fig. le can be explained theoretically, if we 
average the autocorrelation coefficient (258) over all the series, assuming that the ship's 
velocity with respect to the water changes from series to series in 4 random way and that the 


distribution of the ship's velocity obeys the normal law 
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Cc 
wv.) = — em 4] . (261) 


= Jen v CV 


where v is the rms value of the component v,. Using (261) to average (258), we obtain 


+o Co 
Me é 
R(t) = [ exp[- 3 | W(v, av, (262) 
= CO 2 
+o 2 
1 Ab T 1 
= f exp[- 5 + v‘] dv = ; 
on Vv CV a C 
— 00 Cv Co 
l+—tT 
Cc 
a 
or 
Ge ee (263) 
a 
1 + Ba 
db 
where 


If we set T = 0.33 sec, then the function (263) is represented by the curve in Fig. 12, which 
describes the behavior of the experimental points in a satisfactory way. If the scale of 
inhomogeneities is a = 60 cm, the mean velocity of the ship is v~ 128 cm/sec, 1ee. about 


15 km/hour. Unfortunately, no indication of the ship's velocity is given in [23]. 
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THE INFLUENCE OF FLUCTUATIONS ON THE DIFFRACTION 
IMAGE OF A FOCUSING SYSTEM 


INTRODUCTORY REMARKS 


Among the elements of a focusing system there is usually a part which converts a plane 
wave into a spherical wave; for example, this may be a lens or amirror. Since the construc- 
tion of the focusing part is not important, in what follows we shall suppose it to be a lens. 
The field near the focus cannot be calculated using ray theory, since ray theory gives an 
infinite intensity at the focus. By using wave considerations we can find the intensity 
distribution near the focus, the so-called diffraction image. Fluctuations of amplitude and 
phase in the incident Wave give rise to fluctuations of the diffraction image, i.e. the image 
"quivers". When this happens not only do we observe deviations of the intensity from the mean 
distribution, but also the mean distribution itself depends in an essential way on the 
fluctuations in the incident wave. In this regard, two problems arise in the theory of 
focusing systems: 

1) Finding the mean distribution in the diffraction image, 
2) Finding the distribution of fluctuations in the diffraction image. 

Neither of these problems, as far as we know, has been studied before. The problem of 
small fluctuations at the focus of an objective has been studied only in the paper of 
Krasilnikov and Tatarski [10] ; their formulas can be obtained as a special case of our more 
general theory. In this book, we consider both questions without restricting ourselves by the 


* 
requirement that the fluctuations be small . 


a 


The hasic material contained in Part III is given in the author's papers [59,60] . 


a nS SS iS 
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Chapter VII 


GENERAL FORMULAS 


29. The Debye Formula. If the distribution of amplitude A(F) and phase S(F) is known on 


the surface of a sphere behind the lens, then the field near the focus can be calculated by 


Debye's formula [42, DO 


pe sor [ A(F) eo) otk r/F ds. ee 


AF 
Ss 


Here the origin of the coordinate system is chosen at the focus (the center of the sphere). 
The location of the element ds of the spherical surface is given by the vector F, and the 
location of the observation point is given by the vector r. The physical meaning of Eq. (1) 
is clear. In fact, the distance R from the element ds of the sphere to the observation point 


is given by the following formula: 


if the distance r is much less than the focal distance (r << F). The quantity k(F - F-r/F) 
represents the optical path length from the element ds to the observation point. Therefore, 
the contribution of the element ds of the sphere to the field at the observation point is 


equal to 
AP) _tS(F) Zi E-Fer/E) a 


Integrating over the part of the sphere bounded by the diaphragm of the lens, we get a formula 
which differs from (1) only by the factor eae This factor can be explained if we solve the 
problem more rigorously, using the Kirchoff formula. However, dimensional considerations also 
Show the need for the factor Li Xs 


Eq. (1) was obtained by Debye under the assumption that there is no phase aberration (S=0) 


126 


and that the amplitude of the spherical wave along the wave front is constant (A = const). 
The Debye formula remains valid when phase aberration is present and when the amplitude is 
not constant, provided that the change of the phase 5S and the relative change of the amplitude 
A ina distance of a wavelength along the surface of integration is small compared to unity, 
as was shown by Tartakovski [43] . In our case this condition is satisfied. An appreciable 
change in the amplitude or phase fluctuations can occur only in a distance of the order of the 
correlation distance in the medium, and since the wavelength is much smaller than the cor- 
relation distance, the change in a wavelength will be small. 

Using cartesian coordinates with origin at the focus and abscissa along the principal axis 


of the lens, we have 
Fer =xx' + yy' + 2z', (2) 


where x,y,z are the coordinates of the element ds of the sphere, and x',y',Z' are the 
coordinates of the observation point. In this case the amplitude and phase fluctuation on 
the sphere are to be given as functions of the cartesian coordinates, l.e. A(x,y,z) and 


Pie: 


S(x,y,z). The coordinates x,y,z are connected by the relation oy +Z ee: so that only two 


of them are independent variables. Choosing y and z to be independent, we obtain 
Pon 
) 


x = Fo Se Aas? 


In the case of a paraxial lens this relation can be replaced by the approximation 


Pie 
Pep - Z aed se yy) ee (2a) 


In spherical coordinates with the same origin and with the polar axis directed along the 


principal axis, we have 
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Fer = Fr cos y= Fr [cos © cose eae eo cos(P - 6H ‘ 


where F,9,@ are the coordinates of the element ds of the sphere, and r,050, are the coordi- 
nates of the observation point. The amplitude and phase fluctuations on the sphere are to be 
given as functions of the angular variables, i.e. A(0,%) and S(0,@). 

In what follows we confine ourselves to paraxial lenses. Then we can assume that the 
distribution of the amplitude and phase fluctuations on the exit spherical surface of the lens 
coincides with the distribution on the input plane surface, or, in other words, coincides with 
the distribution in the incident wave. Introducing the logarithmic amplitude level by the 


relation A = Ae”, we can write the Debye formula in the following final form 


i -~ikF B+is seus. ee ae 
Dota. *e es [ e exp |ik as, 3 (3) 


FS) 


where A. is the amplitude of the incident (unperturbed) wave. 


50. Statistical Averaging. The fluctuations of amplitude B and phase S are random 
functions of time.. In order to find the average distribution, it is necessary to carry out 
statistical averaging over all possible values of B and 8S. In doing this, we must bear in 


mind that the readings of measuring instruments are usually proportional to Ip |* (to the square 


of the pressure in the acoustical case), so that Ip|* characterizes the mean distribution for 
these instruments. However, we begin by calculating the mean field D: which in the first 
place is simpler, and which in the second place is also needed to calculate the fluctuations 
of the field near the focus. 


It must be emphasized that in order to calculate the average of the factor sore it is not 


enough to know just first order moments (B,S) and second order moments Cghcat we must either 
know the infinite set of moments of all orders or the distribution function of the random 
variables B and S. It can be asserted that the distribution of the random variables B and & 


is normal. In fact, the entire distance which contributes to the values of the quantities B 


and S can be divided up into a large number of segments with lengths of the order of the 
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correlation distance of the refractive index of the medium, and the changes in the quantities 
B and S along different segments will be statistically independent. Therefore, the quantities 
B and S are made up of a large number of independent random variables. Then, according to the 
central limit theorem of probability theory, proved by Lyapunov [ud], quantities like B and S$ 
obey the normal distribution law. 


For the two dependent random variables B and S, the distribution law can be written as 


follows 


f(B,S)dBas = : exp[-(aBo+ 2bBS + eS~)] aBas, (4) 


where 


He Vaesh@ yee Se: (5) 


The parameters a,b,c of the distribution can easily be expressed in terms of the mean squares 
Se; sé and the correlation coefficient Rio’ It follows from (4) that 


"2 


oe eae BS =--—, (6) 
oH oH OH 
whence by (5) we find 
Te oe? ya a) 2 i 
BS See ae or BS Ce eee 


Consequently 
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Ve Vs? Vi - 
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ee. 


Using (6) we obtain 
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The average of the quantity ee can now be calculated easily by using (4): 


+ c 
sere 7% oP H5 (3 5) ads (9) 


ay 2 2 
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We write the double integral in the form 


a 2 a 2 
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The integration over B gives 


rs bas a a 


Then the double integral becomes 
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Since the parameters of the normal distribution satisfy the inequality ec -(o*/ >0O, calcula- 


tion of this integral gives 


Tt le-a-= 2ib 
fol ee 
ac - D 


Finally we obtain 
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It remains to substitute for the parameters a,b and c by using Eq. (8): 


BHsS _ 22, a ee 
e = ox te ee S R,,- (11) 
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Thus, for the mean field we obtain 


— i | -ikF Dey DS ee eu 4) oe a oe int. r/F 
P= yp @ A, exp) $B ds + i B &R,,| |e ds . (12) 
iS) 


However, as remarked, the instruments usually measure Ip|* = p*p. Let us now calculate 


this quantity. By (3) we have 


Am 
O 


p* p = 


ae 


It is easy to show that the factors exp(B, + B.) and exp[i(S, - 85) are statistically indepen- 


5) 
dent. To do so we calculate the average of the product (B,+ B,)(8,- 85); obtaining 


Since any two points in the plane perpendicular to the direction of propagation of the wave 


are equivalent, we have 


and therefore 


(B,+ B,) (S,- Pe) = 0. (14) 


If the fluctuations are independent, the correlation coefficient vanishes. The converse state- 


ment is in general incorrect, i.e., vanishing of the correlation coefficient does not imply 
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statistical independence. However, if the fluctuations obey a normal law (are normally corre- 
lated), then the converse statement is also true. Therefore (14) implies the statistical 


independence of the fluctuations B+ B,, and Si7 S52 despite the fact that the four quantities 


Be ge SBicg hoe, ee. 35 are pairwise dependent. 


cer 
Any functions of the independent fluctuations will also be statistically independent. In 


17 
particular, this applies to the functions exp(B, + B,) and exp[i(S, - S,]. Therefore 


aie = 878, - Se ire) 


(15) 


For the difference in phase fluctuations & = S47 9. we write the distribution law in the fol- 


lowing form 


By (16) we have 


+ oo 
i(S.-S ol 2 7) 
Cr PP neg f ate eggle wee. ae a ees Vs (17) 
V Qa E me e 
~ 0 ot 
Moreover 
2 a =o 
E = (Se 8.) = 8, +8, - 28,8, = 28 (l- R.), (18) 
since gf = gé = a. Therefore 
li 2 
: ee 
ans = = 
- 1 35 - Shay Gib R.) 
‘ (19) 


Similarly, for the sum of the amplitudes 7) = B+ B, the distribution law has the form 


2 
ape ae mut 
W(n)dn = 7. : 3 dn; (20) 
ex n en 
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from which it follows that 


Since 


we have 


—— —s_—2 
3, +B, ; = (1+R, ) 


By (19) and (23), Eq. (13) becomes 
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Eq. (24) can be written in the following way: 


O a 
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(21) 


(22) 


Here the expression Aeen® represents the mean square amplitude of the wave incident on the 


surface of the lens. In fact, the amplitude A on the surface of the lens and the logarithmic 


amplitude level are connected by the relation A = Aen 


tion using the distribution 


Oo 
v(B)aB = —L= exp (- = jas, 
V On Bo CB 
we obtain 
oO Oo ope 


2... 52.28 


so that A = Ae 


Averaging this rela- 


Finally Eq. (25) can be written in the following way: 
— 2) 2 
2 B (R,-1)+8 (R_-1) 
eats A S ; Kk; ke ae 26 
Be ee I] : exp ji R(FL-Fa)+t) dsj de,- (26) 
5 8 


Now it is easy to obtain the formula for the mean square fluctuation of the field near the 


focus. In fact we have 


|Ap|* = |p - o|* = (p*¥ - p*)(p - p) = p*p - p* p, (27) 


whence by (12) and (26) we obtain 


oe =) 7) = 
2 B’(R,-1)4+S (R_-1) (2 .2 oe ee 
Soe [| é “b PT aBere Vases or 2 ds,ds,- (28) 
5S 5 


Ne 


It follows from (28) and (26) that not only the fluctuations but also the mean distribution 
depend on the amplitude and phase fluctuations in the incident wave. 
== 2 
2 2 2B 


It is necessary to make a remark concerning the formula A =A » which implies that the 


of 
energy flow in the incident wave increases with distance. Of course, this result contradicts 
the law of energy conservation, since there are no sources which might increase the energy in 
the incident wave as it propagates. The reason for this contradiction is contained in the 
method of small perturbations (both in the ordinary and modified forms), within the framework 
of which our investigation is carried out. The point is that for the zero order approximation 
we take the plane wave Vee kx, on which the scattered waves are superimposed. Since attenu- 
ation of the plane wave is not taken into account, its energy is regarded as constant, and 
more and more energy from the scattered waves is added to this energy as the path traversed by 
the wave in the inhomogeneous medivm increases. It is this which leads to the energy of the 


resultant wave increasing with distance. However, the scattering actually diminishes the 


energy of the primary wave, i.e., the regular wave field is converted into an irregular field 
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in such a way that the resultant flow remains constant and does not depend on the distance. We 
can make the theory agree with the law of energy conservation if we introduce the normalizing 
e 
-B ; ie : asta : 
factor e and set p'’ =e p in the incident wave. Since the normalizing factor depends 
only on x, it attenuates all the waves emitted by different elements of the lens to the same 


extent and consequently does not affect the distribution behind the lens. Using (11) we 


obtain. 


pt 
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for the average normalized field in the incident wave, from which it follows that the mean 
field goes to zero as the distance increases. The transition to a normalized field is equiva- 
lent to changing the amplitude A to the amplitude A, in Eqs. (26) and (28), which characterize 


the mean distribution and the distribution of fluctuations near the focus. 


31. Various Special Cases. The basic formulas (26) and (28) of the preceding section 
were obtained without special assumptions concerning the form of the correlation coefficient 
for the refractive index and without the restriction that the amplitude and phase fluctuations 
in the incident wave be small. Suppose now that the refractive index correlation coefficient 
has the form N = exp(- Be ae Then it follows from Eqs. (219), (220), (225) and (226) of 


Part II that 


B +S =a; (29) 

- -_ Bg 2 

B(R, - 1) +8*(R, - 1) = afe™’ /* - 1), (50) 
where 

oe EE a. a 


(It is interesting to note that in this calculation the terms containing the integral exponen- 


tial function cancel each other.) The formula (26) for the mean distribution becomes 
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a 2 ff afew 8 72°)-T) oe (4 BG. #7) is, ds. (52) 


If the amplitude and phase fluctuations in the incident wave are small (Bo <<...1; 5 iy the 
exponential factors in (28) can be expanded in series so that 

ne ae = 

B (R_-1)+S (R_-1) 2 2 — — 

e * ‘ - ae oe BR, + g°R . 

Ss 
Then from Eq. (28) we obtain 
2 a2 
LA = = : . 
[o> | ,) ew |t EG,-#,)-2] as,ae, (33) 


If in addition we assume that the refractive index correlation coefficient is N = exp(-r“/a‘), 


then Eq. (43) becomes 


2 
—_—— A @ ole 
rag O - t/a Ue, Ae we soar 
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Chapter VIII 


THE MEAN DISTRIBUTION NEAR THE FOCUS 


Finding the mean distribution near the focus when there are fluctuations in the incident 
wave 1S a more complicated mathematical problem than finding the distribution in the presence 
of regular amplitude and phase aberration. The complication stems from the necessity of cal- 
culating a four-fold integral (two integrations over the surface s) instead of a double inte- 
gral. In order to solve the problem even for limiting cases, it is necessary from the very 
beginning to make the simplifying assumption that the lens is limited by a square diaphragm. 
It seems reasonable to suppose that the field near the focus of the lens will agree in its 
main features with the field of the lens when it is limited by a circular diaphragm. In con- 
Sidering a lens limited by a square diaphragm it is expedient to use a cartesian coordinate 
system with the abcissa along the principal axis of the lens. By (2a) we have 


pa ed 2 2 ey x! ; ; 
= “r= s is = = + = 
Ee ye 5 eee i ee ee 


Since 


2 2 
P= (y, - ¥,)° + (2, - 25)*, 


in cartesian coordinates Eq. (32) is written as 


5 n/2 2 2 
ae “o. rahi exp | a (ol SU tit ~x (35) 
MF a 


-h/2 
ok ne es ee oy 
x exp (3 (y8-veae5-2 BS + (y4-¥5)y" +( Z~Z)Z ; dy, 42, dy 425 


where h is the dimension of the square diaphragm. We shall use Eq. (45) to study the distri- 


Lf 


bution both in the focal plane and along the principal axis of the lens. 


32. The Mean Distribution in the Focal Plane. In the focal plane (x' = 0) Eq. (35) takes 


the form 


2 h/e 2 
4 2 (Veen) tees) 
P*P = po TFT xp =o) 1 <2 x (46) 
XN F a 


-h/2 
ae < 
xX exp (: Ely" (y,-¥p)42" (2,-29)]) dy, dz, dydz,, ; 


The integrals which figure in this expresseion can be calculated in the two limiting cases of 


very small (a << 1) and very large (a >>1) fluctuations in the incident wave. In the first 


case a << l. If we set the first factor behind the integral equal to unity, then Eq. (36) 


simplifies greatly to 


Neat 


a 4 k tings ae te. 21 2 
p"p = S/T fp | WY 'Yyc¥'¥otZ'Z, -2 2) |v, a¥_02, 225. (50) 


-h/2 


The fourfold integral factors into single integrals, which are simply evaluated. Then we 


obtain the familiar formula for the intensity distribution in the focal plane of the lens in 


the absence of fluctuations in the incident wave: 


Z nt Sin 
oe ee 
F 


r) (38) 
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The dimensions Ay' = Az' of the central focal spot are given by the formula 


re ee ees (39) 


lL, =p AS (40) 


The study of large fluctuations in the incident wave can give substantially new results 
only in the case where the dimensions of the lens are large or at least of the same order of 
magnitude as the correlation distance in the medium ( h 2a). For if the dimensions of the 
lens are small compared to the correlation distance (h << a), then the amplitude and phase are 
identical at all points of the lens, and their fluctuations in time, even if large, have no 
effect at all on the diffraction image. Therefore, we consider the case where the two con- 


ditions 


a>»>l (41) 


and 
h2a (42) 


are met. In this case we can use the following approximate representation of the exponential 


funetion: 
2 2 2 2 
Ozig Yn) as (2, - Zr) ‘ ae Yo) + (Z5 Z,) (43) 
SA en eS 2 
a a 
For 
2 2 2 
on 2 ve) + (25 = Z) <a 


two terms in the expansion give a good enough description of the exponential function. For 
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the contribution to the integral (36) made by the exact function or by its approximate expres- 
sion (43) are relatively small (the smaller the larger a). Therefore the form of the function 
for large values of the argument is not important. 


By using (43) Eq. (36) takes the form 


: h/2 
= x a 2 ky' 
PP = 32 [| ox|- crn i) ha Yo) lav,a¥, x (44) 
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We introduce the relative coordinate y = Yi- Yo and the center of mass coordinate 
1 
= oat Yo): In calculating the integral an important role is played by the values of Vy 


which satisfy the condition 


For sufficiently large a 


yf << ner 


Then by (42) 


yo<<h =) 


a fortiori. If this condition is met, we can integrate with respect to y in (44) between the 


limits - ~ and +™. Then we obtain 


h/2 
[| exp|- BN ¥)°+ Ly, - Yo) layay, = (46) 
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The second double integral in (44) is identical with the first. Therefore, we finally obtain 


B¥p = I exp[-p(yt@+ 2'*)], (47) 
where 
2.2 
I= tah ae, 18 
ees : ( 
1 ,Ka 
B = = aa) ; (49) 


Introducing the radial distance if = ie ie in the focal plane, we write (47) as 


p*p ae Te Pr 2 (50) 


se) 
\ 


Thus the intensity falls off monotonically as we go away from the focus. By (49) the effec- 


tive radius of the focal spot rio, = 1/ VB is 


oF 
"eff ~ vo ka ’ (51) 


1l.e., it increases as the fluctuations increase, as is to be expected. The dependence on the 
focal distance and the wave number is the same as in the absence of fluctuations. However, 
the correlation distance takes the place of the dimension of the diaphragm. This is quite 
natural, since in the absence of fluctuations the amplitude and phase are constant over the 
area bounded by the diaphragm, while in the presence of fluctuations, the amplitude and phase 
are approximately constant over an area whose dimensions do not exceed the correlation distance. 

The value of the intensity in the focus is given by the quantity I. As can be seen from 
(48), this quantity decreases as the fluctuations increase. Thus, as the fluctuations 


increase the effective radius of the focal spot increases, while the value of the intensity at 
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the maximum decreases. The spot becomes more "smeared out". However the total flow of energy 
through the focal plane remains equal to the flow of energy through the diaphragm of the lens. 
To confirm this, we calculate the total energy flow through the focal plane. It is propor- 


tional to the integral 


co 
f p*p enrdr. 
O 


Substituting from (50) for p*p, we obtain 


Using (48) and (49) to replace I and B, we finally have 


f p*p enrdr = nmas 
O 


The right hand side of (52) is proportional to the flow of energy through the diaphragm of the 
lens. Thus, the flow of energy through the focal plane is equal to the flow through the dia- 
phragm of the lens. Eq. (47), (48) and (49) satisfy the law of conservation of energy. 

By (39) and (51) the ratio of the effective radius to the dimensions of the focal spot in 


the absence of fluctuations is given by the formula 


, (55) 


For equal intensities in the incident wave, the ratio of the intensities at the focus obtained 


from (40) and (48) respectively is 


I Ta 
i. ae. 54) 
ie ae ( 
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Finally, we consider a numerical example. Suppose that the dimensions of the lens are 
equal to the correlation distance (h = a). Moreover, let Vis? - x, which corresponds (because 
of the statistical independence of phase fluctuations separated by a distance h = a) toa 
phase difference of ex between the extreme rays incident on the lens, or to a one wavelength 
path difference between the extreme rays. Suppose further that the amplitude fluctuations are 
negligibly small compared to the phase fluctuations, i.e., Be << 8°, which is true for the 
region of small values of the wave parameter. Then a@ = 5° + ae ine 10 and the conditions 


(41) and (42) for the applicability of the theory can be regarded as satisfied. Then Eqs. (53) 


and (54) give - 
ie J 
Aye Se gee 


The effective radius of the focal spot coincides with the dimensions of the focal spot in the 
absence of fluctuations. The intensity at the focus is 0.3 times the intensity at the focus 
in the absence of fluctuations. If the amplitude fluctuations cannot be neglected, then tak- 
ing them into account increases the dimensions of the focal spot and decreases the intensity in 
the focus, aS can be seen from Eqs. (53) and (54). Increasing the dimensions of the diaphragm 


leads to the same results. 


33. The Mean Distribution along the Principal Axis. On the principal axis of the lens 


y' = z' = O,and Eq. (35) becomes 


h/2 
Ses Ae (y,-¥5)°+(z,-25)° 
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grals of the form 
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where C(z) and S(z) are the Fresnel integrals. Therefore 


2 eae Gv) Gr (56) 
er ve) 


As we go away from the focus along the x' axis, the intensity decreases in an oscillatory 


fashion and approaches zero. 


Let us consider the case of large fluctuations, when the conditions (41) and (42) are met. 


Using the approximation (43) for the exponential function, we transform (35) into the form 


_ eee 
P*P 2 35 ay; om 5 “a(Iy-¥o) oe 2808-98) favyay, * (56a) 
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Going over to the relative coordinate y= Y47 Vo and the center of mass coordinate 


y= ayy Yo) we can integrate with respect to y between the limits - » and 4 o, in the way 


Shown in the preceding section. The double integrals take the form 


h/2 7 
2 ; t 
f ay, f expl- Sy" +i —3 ray |ay- (57) 
— 0o a. 
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* 
Consider first the evaluation of the inner integral . Introduce the new variable 


Soy 2b, 
where ae ky x" 
n= os = (58) 
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Then the integral becomes 
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(The integration along the line parallel to the real axis can be replaced by integration along 


the real axis.) Now we can easily carry out the integration with respect to y, in Gee 


aE i onn[-(—5) ° of, wy (CE) 


kahx' 
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” See footnote on p. 101. 
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where 


, 2 
Q(z) = : i e* ax 


 (n 


is the probability integral. Eq. (56a) can finally be written as 


ee: ety; 


=y 1 yx! 


where I is defined by Eq. (48), i-e. 


cece 2 kah 
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As the argument z increases, 0(2) increases more slowly than z. 


(59) represents a function which falls off monotonically as we go away from the focus. 
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(59) 


Therefore the distribution 


Chapter IX 


FLUCTUATIONS BEHIND THE LENS 


The theoretical study of the distribution of fluctuations in the diffraction image is an 
even more complicated problem than finding the mean distribution. Even in the case of small 
fluctuations in the incident wave the distribution of fluctuations can be found only in the 
two limiting cases h << a and h >a. However, at the focus itself the fluctuation can be 
calculated for any value of the ratio between h and a. The next two sections are devoted to 


a study of these matters. 


34, The Distribution of Fluctuations. By (34) the distribution of fluctuations in the 


focal plane is given by the formula 


h/e 
a Ao. Cie oa F (2, - Zn) 
Jeol” = 35 [ff feo) SS x (60) 
-h/2 


&k 
x exp]i F Ly'(y4- Yo) + Z (2, - Zn) | dy, dz, dy,,4z,. 


If the dimensions of the diaphragm are small compared to the correlation distance (h <a), 
then the first exponential factor in the integrand can be set equal to unity, and the integral 
(60) reduces to the integral (37). The distribution of fluctuations resembles the mean dis- 


tribution ‘Gione If we take the ratio of the average fluctuation to the average intensity at 


the same point, we get 


2 — ———s 
| 4p | Hs By Be asco (61) 
Ip |* 


The fluctuation of the field in diffraction image depends not only on the amplitude fluctuation 


14-7 


Be but also on the phase fluctuation ge in the incident wave, and moreover in the same way. 
There is nothing surprising in this, since fluctuations of the field p (e.g. the pressure) are 
caused both by amplitude and phase fluctuations in the incident wave. 

If the dimensions of the diaphragm are large compared to the correlation distance (h >a), 


then introducing as usual the relative coordinates y,z and center of mass coordinates A 


we can extend the integrations with respect to y ani z from- » to» . From (60) we obtain 
: n/2 ih 
2 ror yo 
|Ap|“ = dy exp|- +5 + i > ylay x (62) 
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The double integrals in (62) are the same as the integral (46) of the preceding section. The 
difference is only in the coefficient before - in the exponential fies instead of ae 


Therefore we immediately get 


ere 


2 
|4p|“ = Ce , (63) 
where 
2.2 
C= “S a ; (64) 
ak,2 
5 = aa (65) 


The fluctuations fall off monotonically as we go away from the focus, and therefore the dis- 
tribution of fluctuations is quite different from the mean distribution, which for small 
fluctuations is approximately described by the relation (38) and has an oscillatory character. 


By (34) the distribution of small fluctuations along the principal axis is given by the 


formula 
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a 
~h/2 


as 
O 
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In the case h < a the distribution of fluctuations along the axis is the same as the mean 
distribution and is given by Eq. (56), in which it is only necessary to introduce an extra 
factor of a. In the case h >>a, the right hand side of (66) factors into integrals of the 


form (57) and we get a solution like (59), i.e. 


Oo 

|p|" = EC |, 2 (67) 
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The distribution of fluctuations is not the same as the mean distribution, which is approxi- 
mately described by (56). As we go away from the focus, the fluctuations fall off monotonic- 
ally, whereas the intensity in the mean distribution oscillates. 

In the case of large fluctuations in the incident wave (a >>1), we can neglect the term 


<2. 2D = ae : 
exp(-Be 7 go) in Eq. (28) compared to ox |3°(R, 1) +8 Coe vf then Eq. (28) agrees with 
Eq. (26) for the mean distribution. All that was said about the mean distribution for large 


fluctuations in the incident wave remains valid for the fluctuations near the focus. 


35. Fluctuations at the Focus. By (43) small fluctuations at the focus are given by the 


following general formula: 
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from which it follows that the fluctuations at the focus depend both on amplitude fluctuations 


and phase fluctuations in the incident wave. Since at the focus 


In| = 
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we get 
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for the relative fluctuation. Instead of calculating the fluctuation of the field at the 


focus, we can calculate the fluctuations of amplitude and phase. Setting 
; I 
p= a exp[i(p- r+ 9], 


we can rewrite (1) in the following form: 
if 1 is 
ae = oe f Ae ds. 
Ss 
Kquating the real and imaginary parts separately, we obtain 


al: 
a cos 4) = 5 [ A cos §S ds, a sin g = f A sin S ds. 
S Ss 


In the case of small fluctuations these become 
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L 
ad = Vr [ ASds. (72) 
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Averaging the relation (71), we obtain 


—- 1 = ls 
a= Wa | Ads = \F As. (73) 
s 
Denoting the fluctuations of amplitude at the focus and in the incident wave by a' =a-a 


and A' = A- A respectively, we find that 


dl. 
a! = iF f A'ds. (74) 
From (72) we get 
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with an accuracy up to quantities of the first order. Dividing (74) and (75) by (73) we find 


that 


From this we obtain 
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Egs. (76) and (77) were obtained by Krasilnikov and Tatarski [10]. The fluctuation of the 
amplitude of the field at the focus is determined from the fluctuation of the amplitude of 


the incident wave, and the phase fluctuation at the focus is determined from the phase fluc- 


nig ae 


tuation in the incident wave. Only small fluctuations at the focus have this property. The 
amplitude (or phase) fluctuation at any other point of the diffraction image depends both on 
the amplitude fluctuation and on the phase fluctuation in the incident wave. 


By using (30), Bq. (70a ) can be written as follows: 


|p| ‘ we ¥)° ae ee a 
a rrp exp an ae dy, dz, dydZ,- (78) 


tel? -h/2 


This integral can be calculated without making any assumptions about the ratio between the 
dimensions of the diaphragm and the correlation distance. The problem reduces to calculating 


the integral 


where y = Yi" Yo is the relative coordinate. Calculating the inner integral gives 
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Introducing the variable 7 -(n/22) (y°/e), we transform the second integral into the form 


h/a h/a h/2a 
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This integral can be found in tables [39] : 
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Finally, we obtain 
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ie ay Le _ 


for the fluctuations at the focus. The graph of this function is shown in Fig. 13. As the 
dimensions of the diaphragm increase, the fluctuations decrease. This explains the averaging 
effect of a large lens. In particular, this result is supported by the fact that the twink- 
ling of stars is much more noticeable to the naked eye than when observei with a large tele- 


scope. 


Cc 
Oo 
3 


Fig. 143 Dependence of the relative fluctuation of the field 
amplitude at the focus on the size of the diaphragm. 
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Appendix I 


It is apparent from Eqs. (99) and (100) of Part II that the phase fluctuations and the 


amplitude fluctuations are given by the real part and the imaginary part of the expression 


Pe ee 2 
W(x,y,z) = - 2& f os xp |i oP el lent) dbanae (1) 
Ys = EP |E Be Eyl Sagan 


respectively. We shall show that the exact formula (89) actually reduces to Eq. (1) in the 
case of large scale inhomogeneities. 


We rewrite the exact formula as 


ice 
Wes 2) a — ae) 2 exp|ix[e-(x-8)] | CES hdednac, (2) 
where 


cS CREE a Gree weet). 


We assume that the refractive index fluctuations can be expanded in a Fourier integral 


CE apse) ~ ff O(EK,5K,) exp[-i(k7 + K 6 ak dk, (3) 


Then (2) can be rewritten as follows: 


Wi(x,y,2) =- ie Iefpy 2 exp |ixb-(-8)] | P(E, Ky Ks) x (4) 
bY “x 


x exp[-i(k,q + k,b)Bédndtdk dk, = 


19) 


-_ a FEL L exp(ikr) exp (2fka(y-n) A Ks(2-$)]) x 


x exp{-i(K,y pe) = ik(x-e)] B(E,K,»K,)dbdndgdk, dk, . 


We introduce polar coordinates in the 7,€ plane by setting 


Then 


Katy zs n) + K(z - €) = kp cos(e€ - a). 


Integrating with respect to « in (4) gives the Bessel function of order zero: 


x ae 
viz - ike f dé {| exp[-ik(x-&) - i(K,y+k,z)] O(E,K5,K,)dk dK, (5) 
oO - 0 


co 


1 ; 
x f — exp(ikr) I (kp)pdp. 
O 


Using Sommerfeld's formula [45], we expand the spherical wave in cylindrical waves: 


= exp(ikr) : f oo T(X6) ex |- Via (x-)| (6) 
O Vi -k 


(e238 70) 


Then the inner integral in (5) can be written as follows: 
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J o(ke) T(4e) daa. (7) 


F ; ee 


Because of the equality 


f(x) = [ ode f 200 Ij(ke) I (ro) adn, 
O O 


[e@) 


known as the Fourier-Bessel integral (valid in the case where J & 120) [ax exists [e7]), the 


double integral (7) equals : 
2 2 
exp[- VK"- k” (x-€)] 
Koo i 
Finally we write Eq. (5) as follows: 
x + 0 coll 
ai f 7 heey = exp[-i(k ,y + k,2)] x (8) 
O - w 


ij2 
x exp[-ik(x-&) + Hic 7) (x-€) | BCE KK.) Abd dK,» 


Consider now the approximate formula (1). Using (3) to replace p(&,7,6), we obtain 


D x + a 
k ae 
W'(x,y,2) = — f x-€ 7 P(E Kp, K,)dK dK, (9) 
0 - © 
i : r (2 -t)" 
- . Line . 
x f exp |ik ne UW ik jan f exp(ik Diane 6) = ix,¢|dt. 


Af 


We calculate the last two integrals by the saddle point method. To do this, we introduce new 


variables & = E/x, ny = n/x; es y/x. Then the first of these integrals can be written as 


+o e 
— fOy ado Ke, 
x | ex ikx| Te - = 1) Pn: (10) 
ak 
Writing 
2 
f(n,) =i wa _ 2 Gre 
il = Ee) gk 


and bearing in mind that the parameter kx is large (kx => 1). we use the saddle point method 


[27] to evaluate the integral 


+ 00 
i= f exp [kxf(n, )]dn,- ies 


=- © 


We find the saddle point from the condition of /Sn, = Oy. Lees 


O 


J, 7 J K 
1 
is = ce ws 
1 
whence 
Oo __é 
Ny ye as ee a (15) 
Furthermore 
ar 4 
a = 
ony O : a 


Making a series expansion of the function f(€,>n,) near the saddle point, we obtain 


Co 
Co 
O 1/ or 2 2 Ko : 
f =f Sf) Ape ee BEC gs ce ey 
1/y 
1 
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We must choose the path through the saddle point so that 


An = (1 + i)s, 


where s is a real parameter, Then 


C 
a ee mi, a 2 
f(n,) = - aoe Pe) Side ae ae (15) 
al 
Cc + 00 
: eee [ kx =o 

ES afi exp|-ix a - Bae ik ,y] f ex - ee Plas (16) 

: C 
= Bans ex eee —E_)- ik 7 | 
Kx BAe Be eS a es 
Eq- (9) then takes the form 

X + 

Wi(x Vy zZ) =k exp|-i(k ytrR ie a ee + Kk?) (x - E)] x (17) 

anes C 4 Ck. 2 4 

Oo - w 


x BCE Kp Kz )dEGK dK, 


It is now clear that in the case of large scale inhomogeneities (k >> x) the exact formula 


(8) reduces to the approximate formula (17). To see this, it suffices to set 


(62 ey V2 ie in the amplitude and to take only the first two terms of the expansion 
2 2,1/2 1K* 
(k™ ~ kK’) ~ k(1 - = —=) in the phase. 
C “o 


ID? 


Appendix II 


To evaluate certain integrals appearing in the text (Chapters 5 and 6) it is necessary 


to use the identity 


eS 1 


+ o 2 2 
. (y,+ y)*4 (25+ 2) 
silos 27 ale eas 
Pane | 


(Yo7 y,) + (z.- 2) 


1 
x ex) | =} x 
2s, en(S, + S,) 


2 2 
Gres Yo) Pt 2) 


X exP|- ———srey ety) 
2 5+ 55 


which is well known in the theory of probability Rae (The result of convolving two 


Gaussian distributions is again a Gaussian distribution.) First setting 
ul 
2 


in this identity, and then 


N 

Ht 
KO) Vans 

+ 

ron} 


a acta) mes 
tao pe ey 


we can write the following two formulas: 


NO how 
+ 
N 
—” 

nh 


+ oo 
q 2 
1 Ll etal e4 : 
2 * Cia 
h 
+14 a, 2 1 


- © 
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(1) 


(2) 


— iyds = Baa, a) * 
~2ia,, TL oat a 
xX exp]- “+ aus Lt )e 
al - ‘ 
a 
+ 0 
1 im Ce ae (2 + z) 
- e ae cla ‘ 
Lise aA, ~ S a 


1 


x €XP]- dydz = pas ao) x 
ela, Oni a+ 4, 


Going over from exponential functions to trigonometric functicns, we obtain 


7 (2 os y)* + (2 + z)* 
ae C eC 
lise or hae ay 


(A-y)? + (Qeer-2)° (A+ y)> + (242) 


x cos |——————_=—_———_ + sin] ———__—_—_____| x 


ca, 8) 
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(3) 


(D+ y)° + (+ 2)° 
cos 5 x 
“a 
(Q- yh + (B+ et - 2)F| 
+ 2. Sain 
28, 
(A-y)? + (B01 - 2)® 
dydz = 
4, 
cos 
= 2) en(a,- a, ) 
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x COs dydz = x 
2a, ana, + a, ) 
n° + (€ + Lrye ; 1. n + (6 + aa 
x Sin B(a a eee. See 2 re ° 
i.e 2 aa ne a al 


Aiding these equations and equating the real and imaginary parts separately, we find 


= 
(A+ y)© + (24 2)° 
2 ; Pas 
eds TT alo als 
mm 
na, &, 1 
(A - y) + (B44 - 2)? ; 
xX sin dydz Pa ( ea 
c8,, tla, - a, 
xX sin ; = we : sin ie 7 wye 
> 
(a, - 2) an(a, + a) 2(a, + 8, ) 
Cry) +(2+2)° 
= [| cos x 
a 
Lee a. Bn ds 
(A- y)® + (B40 - 2)? : 
ee 2 en x(a, - a) " 
a n(ay 7 
2 es 
x cos : as wre + : cos z aa 
(a, - a.) an(a, + 8.,,) a(a, + a.) 
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Suotracting (5) from (4), we obtain 


X cos dydz = x 
28, ana, - a,,) 
dno + # LN/e 1 ee + Erye 
ae ones "Bala, +a.) | 2fa, +a.) |’ 
aa, - a) ohio ae Le 
: (D+ y)? + (242) 
5 i? sin 3 x 
ig a_a om 
Yas 
(A-y)? + (B+ et - 2)* 
2 ae 
X cos - dydz 
a, ena, + a, 
24 (b 4 4? ; 2 + 47? 
a ee | ee, | Se 
a PS aes * ae “ee 


In the notation of Section 18 of Part II, Eqs. (6) - (9) can be written as follows: 


pp 
[fhe Vero Ea» 
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if ify va + ' Gra] x (22) 
x 5 ue = y)* a (S a ae | dydz = 


5 B(2,- app) + O(a, + <2»0)], 


+ © 
{| ane ie + y)< + (2 7 | x (12) 


ae 2 zs 
« Bylegs (B= wi = G+ 2] ayes = 


as 
a By B, 8.» 9) = §, (a,+ a1) ) 


+ oo 
[| Abe ie 2 y)* + (2 + | x (13) 
x ile ya - ae + (2 + 4° - a dydz = 


al 
= a $,.(a,+ 2.530) 7 ,,(a,- 8.52) 7 


where 


p =n + (t+ Lye. (14) 


In studying the problem of amplitude and phase fluctuations and the correlation between them 


at the receiving point, we must set #£' = 0. In studying the transverse autocorrelation of 


amplitude (or phase) fluctuations at different points we have t' #0. 
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